LONG-TIME BEHAVIOR FOR THE TWO-DIMENSIONAL MOTION OF A DISK 

IN A VISCOUS FLUID 
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Abstract. In this article, we study the long-time behavior of solutions of the two-dimensional fluid- 
rigid disk problem. The motion of the fluid is modeled by the two-dimensional Navier-Stokes equations, 
and the disk moves under the influence of the forces exerted by the viscous fluid. We first derive U'-L'' 
decay estimates for the linearized equations and compute the first term in the asymptotic expansion 
of the solutions of the linearized equations. We then apply these computations to derive time-decay 
estimates for the solutions to the full Navier-Stokes fluid-rigid disk system. 
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1. Introduction 

We consider the system formed by a rigid disk and a viscous fluid filling the whole plane M^. We 
assume that the body initially occupies the disk Bq and rigidly moves so that at time t it occupies a 
domain denoted by B{t) that is isometric to Bq. We denote J-{t) := M? \ B{t) the domain occupied 
by the fluid at time t starting from the initial domain := \ Sq- 
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for t e (0, oo) 


X G T(t), 


(1.1) 


for t G (0,00) 


X G I'(t), 


(1.2) 


for t G (0, 00) 


X G dB(t) 


(1 3) 


for t £ (0, 00) 




(1.4) 


for t G (0, 00) 




(1.5) 


for t G (0, 00) 




(1.6) 


for X £ J^o, 




(1.7) 






(1.8) 



When the fluid has constant viscosity > 0, the equations modeling the dynamics of the system 
fluid-rigid disk read 

du 

— + (n • V)^ - i^Au + Vp = 
div u = 

u{t, x) = h'{t) + uj{t){x - h{t))^ 
hm |u(t, x)| = 

|a;|— >oo 

mh"{t) = - Snds 
JaB{t) 

Juj'{t) = - [ (x - h{t))^ ■ Ends 

JdB{t) 

u\t=o = Uo 

h{0) = ho, h'{0)=£o, u{0)=iOo. 

Here, u = (ui,U2) denotes the velocity-field, p the pressure and T, is the Cauchy stress tensor of the 
fluid: 

J: = -pld + 2uD{u), (1.9) 
where Id is the 2x2 identity matrix and: 

2 V oxe oxk ' 

The constants m and J denote respectively the mass and the inertia of the body while the fluid is 
supposed to be homogeneous, of density 1 to simplify notations. In this work, we assume that the 
solid is homogeneous of density m/vr, implying in particular J = m/2 (we discuss in Section [5] about 
a generalization). When x = {xi,X2) G M^, the vector x^ stands for x"*" = (— X2,xi), n denotes the 
unit normal vector to dB{t) pointing outside the fluid domain, h'{t) is the velocity of the center of 
mass h{t) of the body and a;(t) denotes the angular velocity of the rigid body. Indeed, since B{t) is 
isometric to Bq there exists a rotation matrix 



Se(t) 



cos 61 (t) -sin6'(t)' 
sin6'(t) cos6l(t) 



such that the lagrangian coordinates r/(t, x) associated to the body read: 

r/(t,x) := h{t) + S0(^t){x - ho). 

Furthermore, the angle 6 satisfles 0'{t) = uj{t), and is chosen such that 9{0) = 0. Without loss of 
generality, we assume that Bq is the unit disk centered at the origin: h{0) = 0. 

Given {uo,£q,ujo) G H^(J-q) x x M, satisfying the compatibility condition: 

div no = in Jxi ) ^^o = ^0 + ^jJqx^ on OBq , 

T. Takahashi and M. Tucsnak prove in [22] that there exists a unique global strong solution (n,p, h,uj) 
of ()l.ip - ()1.8p . The construction is based on the change of variable: 

v{t,x) = u{t,x - hit)) , p{t,x) = p{t,x - h{t)) , £{t) = h'{t). (1.10) 
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The new unknowns {v,p) are then defined in the fixed domain [0, oo) x {M.'^\Bq) and system (|l.ip -(fTT 
reads, in terms of {v,p,£,uj): 



dv 
dt 



+ (v 



• Vjv - uAv + Vp = 

div V = 
v{t,x) = £{t) + uj{t)x^ 



lim 



Joj'it) = - 



0, 

Ends, 
x"*" • Ends, 



for {t,x) £ (0,oo) X Tq, 

for (t, x) G (0, oo) X Tq, 
for {t,x) G (0,oo) X dBo, 
for t G (0,oo), 

for t G (0, oo), 



I dBo 

v\t=o = Vo 
m = £o , ^^(0) = wo, 

with 

E = -pld + 2uD{v) 
These solutions verify the following energy decay estimate: 



for t G (0,oo), 
for X G J-Q, 



(1.11) 

(1.12) 
(1.13) 
(1.14) 

(1.15) 

(1.16) 

(1.17) 
(1.18) 



^0 



\v{t,x)\^ dx + {m\iit)f + J\ujit)\'^) 



< 



+ 2v / \D{v){t,x)\'^ drdx 
Jo J To 

VQ{x)\'^dx+ (m|4p + J|tJo| 



To 



Vt > 0. 



:i.l9) 



Relying on this estimate, T. Takahashi and M. Tucsnak prove the existence and uniqueness of a global 
weak solution to ()l.lip - ()1.18p for initial data (t;o,^Oi'^o) such that vq G L?'{Fq) and 



divfo = 0, in F^-, vq ■ n = {l^ + '^ox ) ■ n, on OBq. 
In this article, we aim at studying the long-time behavior of these weak solutions. 



(1.20) 



The long-time behavior of solutions to fluid-structure interaction systems has already been tackled 
in different ways. In a series of papers, several authors study the asymptotics of systems without 
pressure, i.e. where the Navier Stokes equations are replaced by a heat equation [IBl [191 [Ml ES] • In 
this simplified case, the force applied by the fluid on a solid is modeled by the circulation of the normal 
derivative of the velocity-field u on the solid boundaries. In the one-dimensional case in \2A\ I25j. and 
then in several dimensions in I18j. the authors show that the multiplier method introduced in [7] 
to study the asymptotic behavior of solutions to convection-diffusion equations (also applied to the 
porous medium equation in [53]) enables to compute sharp decay estimates and asymptotic expansions 
of solutions up to the second order. Even if the divergence-free condition (|1.12p significantly modifies 
the equations, we will strongly use the results in [TUI [T5]. 

The long-time behavior of solutions for the full Navier Stokes equations in the whole space is also a 
long-standing question that has motivated numerous studies. Applying a Fourier decomposition, M. E. 
Schonbek and M. Wiegner show in |2H [27] that the norm of the Navier-Stokes solution decreases to 
zero, which was a question raised by J. Leray [16]. In [3], A. Carpio obtains a sharp description of the 
pressure, which is given by p = A~^(div (n • Vn)). Representing then the velocity-field by a Duhamel 
formula and using a scaling argument, she computes the development of the solution for long times 
up to the second order. 

Another approach consists in removing the pressure by taking the curl of the momentum equation: 



c^iCurlti -|- u ■ Vcurlti — z^Acurlti = 0, 



:i.2i) 



where curlu is the vorticity of the fluid. Without boundaries, such an equation yields the decay of 
the norms of the vorticity curlu. For cutIuq G L^(M^), such that / curluo ^ 0, T. Gallay and C.E. 
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Wayne prove in [11] that the vorticity behaves as t — t- oo hke the heat kernel 

e 4t 



curl uq 



iirt 

Note that if curluo is compactly supported and integrable, then thanks to the Biot-Savart law, we 
have, for large x, 

f curluo x-^ ( \ 

Of course, this imphes that uq does not belong to L^(M^) UL^(M^) if / curluo 7^ 0. Consequently, this 
theory corresponds to solutions with infinite energy. For instance, in |llj . T. Gallay and C.E. Wayne 
deduce that the velocity behaves asymptotically as t — )• 00 like the Lamb-Oseen vector field: 

J curliiQ x^ 



2tt 

which has infinite energy. 

In a domain with boundaries, system (|1.21|) has to be completed with boundary conditions. When 
Dirichlet boundary conditions are imposed for the velocity-field, one might compute Robin boundary 
conditions for the vorticity but with non-dissipative coefficients. Therefore, working on the vorticity 
seems difficult. In the case of one obstacle surrounded by a viscous fluid (i.e., when Bq is fixed and the 
system reduces to the Navier-Stokes equations in the exterior domain J^q completed with homogeneous 
Dirichlet boundary condition on dJ-'o), the recent works \10\ I12j prove that the first term in the long- 
time behavior of the velocity-field is given by the Lamb-Oseen vector field. Their proofs consist in 
a perturbative argument showing that the decay estimates for the solutions of the Stokes problem, 
which were established in [5l [6l [T7], implies that the nonlinear terms tend faster to zero than the 
Stokes solution. To our knowledge, such decay estimates on the Stokes semigroup are only known for 
fixed domains with homogeneous Dirichlet boundary conditions for the velocity-field. 

The only result considering the long-time behavior of a moving particle inside a Navier Stokes fluid 
is due to E. Feireisl and S. Necasova [8]. However, they assume the whole system to be confined in a 
bounded container and they take into account the influence of gravity. Hence, they obtain completely 
different results with completely different methods. Broadly speaking, they prove that, if the container 
has no vertical wall and contains only one particle, the particle reaches the bottom of the container 
asymptotically in time. 

One of the main steps in |10^ [T2] is to establish — decay estimates for solutions to the linear 
Stokes equations underlying the Navier Stokes system. Such results are known for fixed domains 
with homogeneous Dirichlet boundary condition (see [5l[T7j), but in our case, the linearized Stokes 
fluid-solid system reads: 

dv , „ 

uAv + Vp = 

ot 

div V = 

V = e(t) + uj{t)x^ 

lim \v{t, x)\ = 

ml'{t) = - I Ends 
JdBo 

Jj{t) = - x^ ■ T.nds 
JdBo 

v\t=0 = Vq 

£{0) = io, w(0) = UQ. 

To our knowledge, U'—L'^ (p, q ^ 2) estimates are not available in the literature for solutions of ()1.22p ~ 
()1.29p . To be more precise, in [22] , T. Takahashi and M. Tucsnak construct solutions of ()1.22p - (|1.29p 



for {t,x) G (0,00) 


X To, 


(1.22) 


for (t,x) G (0,00) 


X J"o, 


(1.23) 


for {t,x) G (0,00) 


X dBo, 


(1.24) 


for t G (0,00), 




(1.25) 


for t G (0,00), 




(1.26) 


for t G (0,00), 




(1.27) 


for X G /"o) 




(1.28) 






(1.29) 
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via a semigroup approach. They show that the semigroup is analytic on £^ (to be defined in ()1.30p 
below) in dimension 2. In [2H], the semigroup is also proved to be analytic in the counterpart of the 
spaces 

£<i/5 n /IP in 3D. However, in both papers the subsequent decay estimates are not sufficient for 
our purpose. In the first case, the framework only is considered. In the second case, the authors 
do not obtain sharp decay estimates. 

To state our results precisely we introduce shortly some notations. From a triplet (vQ,io,ujQ) ver- 
ifying ()1.20p . we define a divergence- free vector field denoted Vq on obtained by extending vq by 
Iq + loqx-^ in Bq. Adapted to such Vq, we introduce the functional spaces defined as follows: 

CP = {V € LP{R^), divF = in M^, D{V) = in Bq}, {p G [l,oo]). (1.30) 

When p E [1, oo), we endow these spaces with the norms 



\v\\l,= [ \v\p + - [ \v\p. 



I To J Bo 

It is easy to check that, liV ^ CP ^ then V = ly ^ loyx^ on Bq^ where 

iy = -l V{x)dx, ujv = -l V{x)-x^dx, (1.31) 

J Bo ^ J Bo 

and the normal component of V is continuous across OBq (as in ()1.20p ). In particular, we remark 
that, setting v = V\jrg, there holds: 

Such a space is obviously a Banach space as a closed subspace of LP(M?). A straightforward extension 
of O Theorem III.2.3] yields that CP n C~(]R2) is dense in CP for arbitrary p G (1, oo). For p = 2, the 
space CP is a Hilbert space as the norm is associated with the scalar product: 

f TTi f 

{V,W)c2= V-W + — V-W. (1.32) 

J^o ^ J Bo 

For p G (1, oo) \ {2}, the same bilinear form enables to identify the dual of CP with CP' where p' is the 
conjugate exponent of p. 

Naturally, we endow C°° with the norm: 

For any V G we still have V{x) = iy + wyX"*- in Bq with iy and uoy defined by p.3ip . Hence, 
there holds again: 

\\V\\c^ ~ max{||i;||ioo(jrQ), Iwyl} . 

Our first results concern the Cauchy problem for (jl.22p - (jl.29p in CP and the decay rates of the 
constructed solutions. As in [221 [28], we use a semigroup approach: 

Theorem 1.1. For each q G (l,oo), the Stokes operator of the linear problem ()1.22p - (jl.29p generates 
a semigroup S{t) on C^ which satisfies the following decay estimates: 

• For p G [g,oo], there exists Ki = Ki{p,q) > such that for every Vq G jC." 

\\S{t)VQ\\cP < Ki{vt)v~\\\VQ\\c^ for all t >0. (1.33) 

• If (1 for p G [(/, 2], there exists K2 = K2{p, q) > such that for every Vq G , 

\NS{t)Vo\\Lrir,)<K2{iyt)-"^^-v--^\\VQ\\c. for all t > 0. (1.34) 

• For p G [max{2, q},oo), there exists K3 = K^{p, g) > such that for every Vq G , 



1 1 



\\yS{t)VQU.,,,, < K,{ut)-^^-^~-.\\VQ\\c. for all < t < i, 

°"^^^"^-\ K,{ut)--.\\VQ\\c. for all t > i. 
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Our approach is based on the decomposition of the velocity in spherical harmonics. We show 
that the 0-spherical harmonic verifies a heat equation (without pressure) with dynamical boundary 
conditions. This enables us to compute decay estimates using the multiplier method of Escobedo- 
Zuazua [7] in the same way as in [1H\ I19j . The 1-spherical harmonic is the hardest part. A priori, 
it verifies an equation with pressure and non-standard boundary conditions. However, we show that 
there exists an underlying algebra which enables to reduce this equation to a heat equation (without 
pressure) with dynamical boundary conditions. So, we can again reproduce the method of [7] in the 
spirit of |18l 119] . We do not expand the remainder (i.e., the fc-spherical harmonic for k > 2) m this 
part, as we show it satisfies the Stokes equations with Dirichlet boundary condition on Bq which has 
been studied formerly in several papers O [T7] . 

Going further in the spherical-harmonic decomposition, we are also able to compute an asymptotic 
expansion of the solution to the Stokes system (jl.22p - ()1.29|) for well-localized initial data: 

Theorem 1.2. For all p G [2,oo], and for any Vq G n L^(M^, exp(|xp/4)dx), setting £o = ivo and 



we have 



M = (m- tt)£o, 



]imt'-'/P\\S{t)Vo-Uj^it,.)\\^^^^^^=0, 



t—^oo 



lim t 



M 



'Sit)Vo 



S-KVt 



limsupi < -Foo, 



(1.36) 

(1.37) 
(1.38) 



where 



Uj^{t,x)=V' 



1 



e 4i/t 



27r|x|2 



M-x' 



Before going further, let us emphasize the following important remark, which can be easily deduced 
from explicit computations: provided M ^ Q, for all p > 1, 



< liminfti-i/P||[/_^(t, 



4— S>oo 



\Lp{ro) 



limsupii~^/P||[/^(t,-) 



LP 



< OO. 



The quantity Ai represents the total momentum of the system. Indeed, since any Vq € L^i^) 
satisfying the divergence free condition has mean value. 



Vo -|- m£o 



Vb dx -I- m^o = (m — 7r)£o = -M. 



Therefore, if Ai ^ 0, (|1.36p shows that Uj^ is the first term in the asymptotic expansion of S{t)Vo. 
We deduce also from ()1.37p - ()1.38p that, provided / 0, the solid, whose center of mass corresponds 
to /is(t)vb = fo^S{s)Vo ds, goes logarithmically to infinity and stops turning. 



If = 0, then a careful reading of the proof of Theorem 1 1 . 2 1 vields 



lim sup 



t^oo v|log(t)|l/2 



< +00 



(1.39) 



which implies that the disk converges to a fixed state when considering the linearized equations (jl.22p - 
(|1.29p . Note that the condition A^ = is satisfied in the following two cases: 

• m = vr, that is the case of a solid having exactly the same density as the fluid. 

• ^0 = 0, that is the case of a solid whose center of mass has zero initial velocity. 

Thus, when A^ = 0, we expect a different behavior as t — )• oo of the solutions of the Stokes system 
(fT:22]l - (fL29]l . 
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With the decay estimates of Theorem II .11 at hand, we study the long-time behavior of solutions to 
the Navier Stokes system (|l.lip - (|1.18p . We prove that, for small initial data, such solutions satisfy 
decay estimates similar to the one of the solutions to the Stokes equations: 

Theorem 1.3. Let q G (1, 2]. Then there exists Xo{q) > such that, for all initial data Vq € H £^ 
satisfying the smallness assumption 

\\Vo\\c^<Hq), (1-40) 
the unique weak solution V of (jl.lip - (jl.lSp with initial data Vq satisfies the following decay estimates: 

• for all p G [2,oo), there exists H{p,q,V()) such that 

snp{t-^'-v\\V{t)\\cp} < H{p,q,Vo). (1.41) 

t>o 

• there exists H£{q,V()) such that 

suv{t-^\£vm < He{q,Vo). (1.42) 
t>o 

Besides, the function q i— t- Xo{q) can be chosen as an increasing function of q £ (1,2] which goes to 
zero as g — )• 1. 

The proof of Theorem ll. 31 consists of two steps. First, we consider the case q = 2. Following the idea 
developed by Kato in [T3] , we construct successive approximations Yn which verify the decay estimates 
(|1.4ip uniformly in n for p = 2, p = 8 and ()1.42p . Next, we pass to the limit to get a solution to the 
Navier-Stokes equations with such a time-decay. To reach p € [2, oo) we use a bootstrap argument 
based on the Duhamel formula as in [3]. We then develop the case q G (1, 2) by showing that estimates 
(|1.4ip are satisfied uniformly by the sequence Yn- 

Using then a bootstrap argument allows us to quantify the proximity of the solution of the non-linear 
system (frTT]) - (fLT8]) and of the hnear system (fr22]) - (fT:29]l : 

Theorem 1.4. Let q G (1,2]. Taking Xo{q) > as in Theorem \L3l for any Vq G R £^ verifying 
(|1.40p . the unique global solution V of ()l.lip - (|1.18p with initial data Vq verifies: for all p G [2,oo) 
there exist constants C{p,q,Vo) > for which: 

snp{t'~-p\\V{t)-S{t)Vo\\cp} <C{p,q,Vo), i/gG (1,4/3), (1.43) 



t>2 



t'-l 



s^P<^T7^\\yit)-Smo\\cp} <C{p,q,Vo), ifq = 4/3, (1.44) 



t>2 ^log(t) 

2 11 



suv{t-^-^~p\\V{t)-S{t)Vo\\cr'}<C{p,q,Vo) if qe {A/3,2]. (1.45) 

t>2 

Similarly, there exist constants Ci{q,Vo) > such that 

snpt^^Wvit) - is{t)Vo\} < Ce{q,Vo), if q £ (1,4/3), (1.46) 

sup\j::^\ivit)-esit)Vo\}<Ce{q,Vo), ifq = A/3, (1.47) 



t>2 



.log(i)' 



sup{tt"^|^y(t)-%)yJ}<Q(g,yo). i/gG (4/3,2]. (1.48) 

t>2 

Let us comment the fact that if the initial data Vq belongs to n for some q G (1,2) and 
satisfies the smallness condition (jl.4Up . the U'-noira of the difference between the solution of the 
complete non-linear system (|l.lip - (jl.l8p and the linear one, given by S'(t)Vo, decays faster than the 
a priori decay estimates predicted by Theorem ll.il Indeed, we check easily that 1 — | > \ ~\ fo'^ ^'^Y 
g > 1 and that | - i - i > ^ - | for any g < 2. 

Combining Theorem 11.41 and Theorem 1 1 . 1 1 and taking Aq = Ao(5/4), we can guarantee that, for all 
q G (1, 2] and ah Vq G n £^ satisfying ||Vb||£2 < Aq, we have 

1 _i 

supft? p ||l^(t)||£p} < OO. 

i>2 
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Indeed, for q > 5/4, it relies on Theorem 11.31 and the fact that q i— t- Ao(g) is increasing. For q G 
(1,5/4], it is a simple combination of Theorem 11.41 (for ||Vo||£2 < Ao(5/4) because Vq G with 
the decay estimates of Theorem 11.11 (because Vq ^ C^). If Vq satisfies the further assumption Vq £ 
n -L^(M^, exp(|j;p/4) drc) with ||Vo||£2 < Ao(5/4), we can also combine Theorem 11.41 with Theorem 
11.21 yielding that, for all p > 2: 

SViVt^~h\V{t)\\cv < oo. 
t>2 

In all these cases, we obtain that the solution to the Navier Stokes system thus decays with time (at 
least) as fast as the solution to the Stokes system. 

The paper is organized as follows. In next section, we collect some preliminary results. We explain 
the decomposition of the velocity-field in spherical harmonics. We then compute the different equations 
satisfied by the different modes of the velocity-field and we end up the section by several elliptic lemmas 
that will be used further. Section [3] is devoted to the proof of Theorem 1 1 . 1 1 and Theorem 11.21 Section 
m contains the proof of Theorem 1 1 . 3 1 and Theorem 11.41 The article ends by some comments and open 
problems, in particular the lack of the first asymptotic term of the Navier-Stokes solution. 

Notations. In the whole article, we use classical notations for function spaces. The symbol 
LPiU^dfi) stands for the Lebesgue space with respect to measure d/x defined on an open set Q, C M". 
If d/x is the Lebesgue measure, we drop d^u. Sobolev spaces are denoted by H'"^{Q), m £ Z. Further 
notations for function spaces are introduced along the paper. We shall use extensively symbol L 
in different fonts (such as C^, This will correspond to variants of Lebesgue spaces. The only 

exception concerns ^d^) (resp. -2'c(^ which represent the Banach space of continuous linear 

operators from a Banach space X to itself (resp. a Banach space X to another Banach space Y). 

In what follows, we will use capital letters to denote functions defined on M?, as we did for the 
velocity V above, and denote by the corresponding small characters the restriction on Fq. To be more 
precise, for V, W, Z (• • • ) defined on M?, functions v, w, z denote the corresponding restrictions of 
V, W, Z on Fq and iy, (-wAz the mean value of V, VF, Z on Bq. In the sequel, when considering 
functions W, Z {■ ■ ■) which are constant on Bq, we will identify them with the couple (wjiw), {^Az) 
(• • • ) and write W = {w,(.w)^ ^ = {z,iz)- In the case of the velocity V in C^, the restriction of V is 
iv + i^yx"*- and thus we also identify V with the triplet {v,iv,(^v) and note V = {v,£v,ujv)- 

2. Preliminary results 

We first recall how the Cauchy problem for ()1.22p - ()1.29p has been tackled in [22]. Formal energy 
estimates imply that, for a sufficiently smooth and localized initial data, V{t) G £^ for all t. In 
this framework, system ()1.22p - ()1.29p reduces to the abstract ODE: dtV + AV = 0, where A is the 
unbounded operator with domain: 

V{A) = {V e H\R^), G H'^iTo), divy = in M^, D{V) = in Bq} (2.1) 

such that for V eV{A), AV := F2AV with 

-uAV in To , 

— / D{V)nds+'^x^ [ ■ D{V)nds{y) in , 

and where P2 is the orthogonal projector from L^(M^) onto C^. 

For arbitrary p G [l,oo), is a closed subspace of L^(M^). Hence we can define Pp the projector 
operators from ^^(M^) onto £p which coincide with P2 on ^^(M^) n ^^(M^) (see e.g. ^). These 
projectors are obviously continuous and satisfy WpV = ¥qV for all V £ D L"^. In what follows, we 
omit the index p. We emphasize that the pressure does not appear in the abstract ODE. But, once 
a solution V is constructed, one shows the existence of a pressure such that ()1.22p - ()1.29p holds true 
(see the proof of Corollary 4.3 in [22]). According to this, for sake of simplicity we omit to mention 
the pressure when considering solutions of ()1.22p - ()1.29p . 
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Proposition 4.2 in [22] shows that A is a self-adjoint maximal monotone operator. Therefore, 
applying Hille-Yosida theorem (see e.g. [21 Theorem 7.7]) yields global solutions to (jl.22p - (|1.29|) for 
arbitrary initial data Vq £ C'^. Furthermore, there holds: 

\\V{t)\y<\\Vo\\c2, VtGM+, 
and there exists a constant C such that 

\\dtvmc^ = \\^vimc^<j\mc^. 

Using the identity 

^\\D{V)\\l,^^^^ = {AV,V)c2 
for V G T>{A), see [221 P-61], Lemma 4.1 in [22] implies 

l|VV(t)|li2(^„)<-^||yo|lz:2. 

Hence, previous results in [22] imply Theorem 11.11 when q = p = 2. 

To generalize this result to arbitrary values for p and q, we provide here an original decomposition 
oiV{t). 



2.1. Spherical- harmonic decomposition of spaces. To motivate the spherical-harmonic de- 
composition of C^, assume for instance that V is smooth and denote {iv,ujv) G x M the only 
pair such that V{x) = ly + wyx"*" in Bq. As V is divergence-free, there exists ^ G C°°{M?) such that 
V = V-^^. Fixing §(0) = yields: 

^(^) = 2'^^l^l^ + Iv ' x'^ , in Bq . 
Consequently, introducing radial coordinates (r, 9) and expanding \I' in Fourier series: 

oo oo 

§(r, e)=Y^ ^fc(r) cos(A;^) + ^ $fc(r) sin(fc^), V (r, 9) G (0, oo) x {-tt, tt) , 
k=0 k=l 

we observe that, setting iv,i = • ei, iv,2 = iy ' ^2- where ei and 62 is the canonical orthonormal 
basis of M^: 

^o(r) = ^r2, ^^f-^(r) = ev,2r , $i(r) = -V,ir, VrG(0,l), 

^k{r)=0, ^>fc(r) = 0, VA;>2, VrG(0,l). 

Then, informations on uy and iy contained in the zero and first modes of ^ respectively, so that 
these modes are handled separately from the others. In particular, we focus on dr^o, $1, ^'i, that we 
denote by M^, <I>,^ respectively and regroup the other terms into a remainder. In what follows, we 
still denote {r,9) radial coordinates and introduce {er,eg) the associated local basis. Accordingly, we 
denote by Vr and Vg the radial and tangential components of a vector V. To state our result precisely, 
we also introduce, for p £ [l,oo], the set 

= {V eCP,V = 0on Bo}. 

Though this space contains functions defined on M^, we will often identify the elements of La{J'o) with 
their restrictions on J^q. 

Proposition 2.1. Let p G [l,oo] and V G C^, then there exists a unique {W,"^ ,^,Vr) such 

that: 

(i) V{r,9) = W{r)min{r,l)ee{9) + V^['^{r)cos{9)] + V^[^{r)sin{9)] + VR{r,9), 

(ii) W = W{r) G LP((0, 00), rdr) , and W is constant on (0, 1); W{r) = = uy for r G (0, 1). 
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(iii) ($,^>) = (^(r),^>(r)) G VF/;f(0,oo) with 

p 



f'OO 






p 








+ 


Jo 




r 




r 



|a,.^(r)|2 + |a^$(r)|^ 



rdr < oo , 



and the functions ^ /r,^/r,dr^ ,dr^ are constant on (0,1); ^{r)/r = dr^{r) = £2 = ^v,2 and 
$(r)/r = a,.^(r) = -ii = -£v,i for r G (0, 1). 
(iv) Vfj = Vfj{x) G La{J-o) and the following identities hold true: 



2tt 



VR{r,e)-er cos{e) dB 



2-K 



VR{r,e)-erS{n{e) d9 



2tt 



VRir,e)-ee 



0, Vr G (1,00) . (2.2) 



Furthermore, there exists a constant C (p) depending only on p such that 

ll^llLP((0,oo),rdr) + ll^ij||LP(M2) 

+ l|c^r*||LP((0,oo),rdr) + ll^/''llLP((0,oo),rdr) 

+ l|c^r^||LP((0,oo),rdr-) + ll*^/'^llLP((0,oo),rdr) < C'(p)ll^lkp • 

There exists also a constant C (p) depending only on p so that conversely: 

\\V\\cP < C'(p)(^||T^||ip((o,oo),rdr) + 1 1 ^i?J I LP (M2 ) + || 9,.* ||LP((o,oo),r-dr) + ll^/''llLP((0,oo),r-dr) 



(2.3) 



+ l|c^r<^*||LP((0,oo),rdr-) + II ^/^^ II LP((0,oo),rdr)) 



and 



(2.4) 



\^y\\LP{To) < C'(p)(^l|9rT^||LP((l,oo),r-dr) + II W^AIlLP((l,oo),rdr) + II VVrJ|lp(Jo) 

+ \\drr^\\LP{{l,oo),rdr) + l|f^r^/?'||LP((l,oo),r-dr) + 11^/?^^ llLP((l,oo),rdr-) 

+ \\drr^\\LP{il,oo),rdr) + W^r^ /r\\LPi{l,oo),rdr) + II '^'/''^ II LP((l,oo),rdr) ) • (2-5) 



Proof. Let p G [l,oo]. We first note that, given Vr G Lct(-Fo), it is possible to define by duality the 
functions: 



r I— )■ 



/ VR{r, 9) ■erCos{9)de, r^ [ VR{r,e) ■ ersin{9) d9 , r^ [ VR{r,9)-eed9 
Jo Jo Jo 



on (0, oo). This yields Lj^^{l,oo) functions which might satisfy (|2.2p . Also, once <I>,^' and VF, Vr are 
constructed with the regularity of (ii)-(iv), then (i) yields: 



Vr{r,9) = ^^sm{9) ^-^ cos{9) + VR{r,9) ■ Cr , (2.6) 

Vg{r,9) = W{r)mm{l,r) + dr'^{r)cos{9) + dr<^ir)sm{9) + VR{r,9)-ee. (2.7) 
This implies (g^D and (l23D . 

To prove existence and uniqueness of W,^,^, we assume V G C. With this further assumption. 



identities (|2.6p and ()2.7p together with ()2.2p imply that the only possible candidates W, ^' are the 
following functions: 



W{r) :-- 



$(r) := 



1 



27r min(l, r) 

r27r 



2n 



Ve{r,9)d9, 



Vr{r,9)cos{9)d9, 



^(r) 



27r 



Vrir,9)sm{9)d9 , 



(2.8) 
(2.9) 
(2.10) 



Differentiating the formulas (|2.9p - (j2.10p and recalling that V is divergence-free then yields: 

1 1 /"^^ 

d^^{r) = - Ve{r,9)s\ii{9)d9, dr^{r) = - Ve{r,9)cos{9)d9 , 

Jo ^ Jo 
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(where these identities have to be understood in the sense of T>'{0, oo)) and we have then 

Vn := y- (VFmin(l,r)e0 + V^[^cos(0)] + V^[$sin(0)]). (2.11) 



In the ball, we deduce from these definitions and from V = iy ~^ toyx that for all r E (0, 1): 

W{r)=u;v, ^{r)/r = dr^ir)=iv,2, <^ir)/r = OMr) = -iv,i, Vr = 0. 
From the definition ()2.8p of W, Jensen inequality implies that: 

\W{r)\P<—[ \Ve\P{r,e)de Vr > 1. 
27r Jo 

Combining with the remark that W{r) = ujy for r < 1, we obtain there exists a constant C for which 

WW 

Lp((o,oo),rdr) ^ C'(p) || F || £p . Similarly, we prove that 

l|f^r*||LP((0,oo),rdr) + ll^/''llLP((0,oo),rdr) < ) 
l|5r$||LP((0,oo),rdr) + II ^ A II LP({0,oo),rdr) < C{p)\\y\\cP ■ 

Finally, straightforward computations yield that Vr is divergence-free, vanishes in Bq and satisfies 
(|2.2p . As, combining previous estimates and ()2.1ip also yields that Vr E we conclude that 
Vr E La{J-o) and that ()2.3p holds true. This ends the proof of Proposition 12.11 □ 



Of course, Proposition 12.11 and Theorem 11.11 hold true if we replace the norms || • ||£p by || • ||lp(]r2)- 
We have chosen to keep the notations of Takahashi and Tucsnak [22], where they prove that A is a 
self-adjoint maximal monotone operator for the scalar product (|1.32p . 

Let us also emphasize that in Proposition 12.11 all the functions (W, dr^,'^/r,dr^,^/r, vr) defined 
on are constant on Bq, so that we can identify these extensions with the pairs given by their 
restriction to J-q, denoted will small caps, and their mean value on the ball Bq, denoted i. For 
instance, we will write W = {w,iw)- Moreover, in all the text, we will identify iwit) = ^w{t)- 



2.2. Decomposition in spherical harmonics of the Stokes semigroup. In the rest of this section 
and in Section^ we only consider smooth initial data, namely Vq E £^nC^(M^). Indeed, it is sufficient 
to show Theorem 11.11 for smooth initial data, because H C^(M^) is dense in C^, for the norm 
with q E (1, oo). So the estimates (|1.33p - (jl.35p could be extended, thanks to the linearity of the Stokes 
system. 

In this paragraph, we prove that the spherical-harmonic decomposition of is well-adapted to 
compute solutions of (jl.22p - (jl.29p . We prove: 

Proposition 2.2. Given Vq E £^nC^(M^), the spherical-harmonic decomposition provided by Propo- 
sition [2n\ of the unique solution V E C([0,oo);£2) of ([122]) -((r29|) satisfies: 

• W = {w,iw), w/iere w E C([0,oo),L2((l,oo),rdr)) nC°°((0,oo) X [l,oo)) verifies: 

dtw + V (^-^drirdrw) + ^u?j = /or (t, r) E (0, oo) x (1, oo); (2.12) 
w{t,l)=iwit) /ortE (0,oo); (2.13) 
i'wit) = ^idrwit,l)-wit,l)) /ortE (0,oo); (2.14) 

• = (drip^h) and ^/r = (^/r,^2), where drip,ip/r E C([0, oo); L2((1, oo), rdr)), ip E 
C°°((0,oo) X [l,oo)) and there exists a pressure qi E C°°((0,oo) x [l,oo)) satisfying drqi E 
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C((0, oo); oo), rdr)) such that: 

dti^ + u ^-^dr{rdrip) + ^ip^ = -rdrQi /or (t,r) G (0,00) X (1,00); (2.15) 

dtdri^ + udr (^-^drirdri^) + ^i^^ = /or (t, r) G (0, oo) X (1 , oo) ; (2.16) 

V'(t,l) =5^V(i,l) =^2(t) /or tG(0,oo); (2.17) 
m /I 1 \ 

-i2{t) = -qi{t,l)-i^{—dr{rdri^) + ^^){t,l) /or tG(0,oo); (2.18) 
vr \ r r^ y 

• dr^ = (SrV', — ^i) a^c^ ^/"T = {(p/r,—£i), where {dr(p,f/r) G C([0, oo); L^((l, oo), rdr)), 99 G 
C°°((0,oo) X [1,00)) and there exists a pressure pi G C°°((0,oo) x [l,oo)) satisfying drPi G 
C((0, 00); L^((l, 00), rdr)) such that: 

dt(p + u (^-^dr{rdr(f) + ^(f^ = rdrPi /or (t, r) G (0, 00) x (1, 00); (2.19) 
dtdrif + i^dr ( — -dr{rdr(p) + \(pj = — /or (t, r) G (0, 00) X (1, 00); (2.20) 



r 

^(t,l) =a,^(t,l) = -4(t) /ortG(0,oo); (2.21) 
m /I 1 \ 

-£[{t) = -pi{t,l) + i^{—dr{rdrip) + ^ip){t,l) /or t G (0,cx)); (2.22) 
vr y r r^ y 

• Vr = (-URjO), wherevR G C([0, 00), L2(J"o))nC°°((0, 00) x Jb) and there existspR G C°°((0,oo)x 

- uAvR + VpR = for {t, x) G (0, 00) x ; (2.23) 
div VR = for {t, x) G (0, 00) xTq; (2.24) 
?;jj(t, x) = for {t, x) G (0, 00) x OBq. (2.25) 

We postpone the proof of Proposition 12.21 to Appendix |X1 This mainly consists of tedious compu- 
tations. 

The main interest of Proposition 12.21 is that it reduces the study of the Stokes semigroup to the 
study of scalar equations for the modes involving non-trivial boundary conditions and one Stokes 
equation with homogeneous boundary conditions on the obstacle. Indeed: 

• System ()2.12p - ()2.14p is a scalar heat equation with dynamic boundary condition. 

• System ()2.23p - ()2.25p is a Stokes equation with a fixed obstacle and Dirichlet boundary condi- 
tion. 

• Systems ()2.15p - ()2.18p and ()2.19p - ()2.22p are similar one to each other. Actually, {ip,pi,£i) 
solves (I239|) - (l2:22| ) if and only if (-99,^1,^1) solves (l2T^ - (l2T8D . System (l2T5|) - (l2T8] ) in- 
volves two scalar heat equations ()2.15p - ()2.16p which contain the term qi reminiscent from the 
pressure. It also involves intricate boundary conditions ()2.17p - ()2.18p which couples Dirichlet 
{il^{t, 1)), Neumann {drip{t, 1)) and dynamic (see ()2.18p ) boundary conditions. 

Whereas systems (|2.12p - (j2.14p and (j2.23p - (|2.25p are classical and widely studied in the literature, 
systems (|2.15p - (|2.18p and (|2.19p - (|2.22p do not seem known and are the main challenge of our study. 
Actually, we show that systems (|2.15p - (|2.18p and (|2.19p - (|2.22p reduce to a heat equation with dynamic 
boundary conditions. Concerning ijj for instance, our strategy consists of removing the pressure term 
and reduce ()2.15p - ()2.18p to a scalar equation for the new unknown: 

Z{r) := dr^{r)+'^^ = -dr[r^{r)\, VrG(0,oo), (2.26) 
r r 

which, in particular, is a constant function on the ball Bq, denoted by and for which we have 

Z{r)=lz VrG(0,l), £z = 2£2. (2.27) 
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Note that, using the definition ()2.26p of Z and the fact the ^'(r)/r = I2 on the unit ball, identity ()2.27p 
immediately implies = ^2 on the unit ball, thus being completely compatible with the boundary 
conditions pTT]) . 

Indeed, using this new unknown, we get: 

Proposition 2.3. Given Vq G £^nC^(M^), let {W^ <I>, Vr) he the spherical-harmonic decomposition 
given by Proposition \2.1\ of the solution V G C([0,oo);£^) of (|1.22|) - (|1.29|) . Then, setting Z = {z,£z) 
as in (or Z = -dr{r^)/r), 

• z G C([0, 00); L2((i, 00), rdr)) n C°°((0, 00) x [1, 00)) 

• (z, £z) is a solution to: 

dtz-v(drr + ^dr^z = Q /or (t, r) G (0, oo) X (1, oo); (2.28) 

z{t,l)=lz{t) forte {0,00); (2.29) 
4(t) = aQvdrz{t, 1) for t G (0,oo); (2.30) 

with 

ao = — . (2.31) 

TT + m 

Proof. Up to a change of sign, we focus on Z = dr{r^)/r. Thanks to the regularity proved in Propo- 
sition [2?2l we have {drtp,ip/r) G C([0, cxd); L^((l, 00), rdr)). Consequently z = drip + tp/r enjoys the 
same regularity. The smoothness of z is straightforward. 

Differentiating ()2.15p with respect to r and subtracting ()2.16p . the pressure qi{t) satisfies, for each 
time t > 0: 

dr{rdrqi) + ^gi = 0, for r G (1, 00). 

Hence qi{t,r) = ai 

(t)r + M). Of course, the condition drqi G L^((l, 00), r dr) implies that: 

gi(t,r) = ^, (2.32) 



and therefore, for all f > and r > 1, 

qi 

r 

With this identity, the pressure can be removed simply by adding (|2.16p to 1/r times (|2.15p : 



-drqi 



dt 

Using (I22S1), 



dr + -]^lj 

r 



+ '^{dr + —] (--drirdri^) + 1^ ) = 0. 



r 



drZ = - { --drirdrip) + ] , (2.33) 



and the new variable z in ()2.26p solves (|2.28p . 

Concerning the boundary conditions, ()2.17p reads as 

Z{t,l) = 2i2{t)=iz{t) 

and, using ([2:33]) and (l2:32]l . (l2T8]l yields 

777 

-£'2it) = -Pi{t) + udrZ{t,l). 
TT 

Moreover, still using (f233]l and ([232]), (|2T5]l for r = 1 and ([ZTTP gives 

e'2{t)-udrZ{t,l)=f3l{t). 

Combining the previous equations, {z,£z) solves (|2.28p - (|2.30p . □ 
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Remark 2.4. In what follows, given V = {v,i,ui) a solution to ()1.22p - ()1.29p on (0,oo), we keep the 
convention: 



Zq,{t,r) :=a,^'(t,r) + ^^, Z^{t,r) : = 



, V(t,r) G [0,00) X (0,00) . 



We emphasize that, for t > 0, V{t, ■) has continuous normal and tangential traces through dB^, and 
thus then Z^(t, •) and Zi^[t, •) have continuous traces through the interface r = 1. 

Remark 2.5. As we recalled in the introduction, the classical approach would rather consist in the 
elimination of the pressure in the Navier Stokes system by taking the curl of the Navier Stokes equation, 
yielding that way an equation for the vorticity of the velocity-field. But this is not the method we 
choose here. Indeed, in an exterior domain, one should complete the vorticity equation, and this would 
yield non-dissipative boundary conditions of Robin type. 

2.3. Some elliptic problems. To conclude this section, we prove some technical lemmas that will 
be useful later on. Indeed, in order to compute the decay of the Stokes semigroup, we study the decay 
of solutions to the heat equation ()2.28p - (|2.30p . This gives the decay of the new unknown z whether 
it is computed with respect to ip or ^. However, to our purpose, we need then to invert the definition 
of z in order to get also the decay of (p and ^ in suitable spaces. This is the content of the following 
proposition: 

Proposition 2.6. Given p G (l,oo] and {z,i) £ L^((l, 00), rdr) x M, there exists a unique ij) € 
WjJ^(l,oo) solution to the following boundary value problem: 



drtpir) + 



r 

m 



z{r) 



for r G (1, 00) 



and there exists a constant C (p) depending only on p for which: 



\9rtp\\LP{{l,oo),rdr) + 



LP {{1, 00), rdr) 



< CiP) (ll^lll,f((l,oo),rdr) + 1^1) 



(2.34) 
(2.35) 

(2.36) 



Proof. Let p £ (l,oo] and {z,£) satisfy the assumptions of Proposition 12.61 It is straightforward that 
the unique solution to (|2.34p - (|2.35p reads: 



^(r) 



i 1 
- + - 

r r 



sz(s) ds , Vr > 1 



If j3 = 00, we establish easily (j2.36p from this formula. If p S (1, 00), up to a regularizing argument 
we skip for conciseness, we multiply (I23i]l by \^\P~^i;/rP-^ on for arbitrary R > 1: 



R 



z|^| 



P-2. 



rP- 



■r dr 



> 



\mr)\p] 


R 




j.p-2 


1 




1 




p 


R 


p 


j.p-2 




1 


_\ 


p\p 








^ + 2 






p 







{p-l) 



dr 



1 pr 



p-2 



dr 



R 



rP 



-r dr 



211-i 

P 

R I, 



R 



rP- 



■ dr . 



P 



rP- 



dr . 



Hence, for all p G (1, 00) 

p 



LP{{l,R),rdr) 



< C'(j') ||2|lLP((i,oo),rdr) 

< C{p) ||^|lLP((l,oo),rdr) 



p-l 



+cipm^ 



LP{{l,R),rdr) 



LP'((l,R),rdr) 

+cipmp. 
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This yields 



< C{p) (|k||LP{(l,oo),rdr) + 



LP{(l,iJ),rdr) 

Letting then — )• oo we obtain (|2.36p . 

Let us now state another ehiptic estimate that wih be useful in the following: 



□ 



Proposition 2.7. LetpG (l,oo)\{2} and assume that z G LP{{l,oo),rdr) and drZ £ L'''{(l,oo),rdr). 
There exists a constant C {p) depending only on p such that: 



< C(p) (||5rZ||ip((i^oo),rdr) +£p\z{l)\) 



(2.37) 



LP((l,oo),rdr) 

where Sp = 1 if p > 2 and Sp = if p < 2. 

Proof. As z belongs to W^^^{l,oo), we infer that it is continuous and we integrate by parts on 



1 



LP{(l,R),rdr) 



< 



p-2 
1 



R 



1 



z{l) 



\j.p-2 

z{Rr 



dr 



1 


\ W 1 


R p 


\z\P"^zdrZ 


p-2 


j.p—2 


1 %-2 j 


1 rP-1 



V 



p-2\''''" RP-^ 
Then, the following depends on the sign of p — 2: 
• if p > 2, we directly have that 
p 



\p-2\ 



\\drZ\ 



rdr 



LP((l,oo),rdr) 



p-1 

LP{{l,R),rdr)' 



LP{{l,R),rdr) p - 2 



< -^\Z{1)\P + ^^WdrZ 



p-2' 



rZ\\LP{(l,oo),rdr) 



P-1 

LP((l,R),rdr)' 



which gives (|2.37p with Sp = 1. 
if p < 2, we get 



p-i 

LP{{l,R),rdr)' 



1 \z{R)\P p 

LPiil,R),rdr)- 2-p RP~^ ^ JZT^W^^r z\\ LP ((1,00), rdr) 

To establish (j2.37p with Sp = 0, it is sufficient to find a sequence Rn — )• oo such that 
{\z(Rn)\PRn) tends to zero. This can obviously be done since r i— )• r|z(r)|P is assumed to 
belong to L-'^(l,oo). 

□ 



We finally provide elliptic estimates that will be useful when getting estimates on the 0-mode: 
Proposition 2.8. Let p G (l,oo) and assume that w S LP {{l,oo), rdr) satisfies 

for r E (1,00) ; 



drw{r) w{r) 
drrw{r) + ^ = /(r) 



drw{l) — w{l) = a, = ^1 

for some f G -L^((l, 00), rdr), a, b in M. Then, there exists a constant C{p) depending only on p for 
which: 



|9r-r-^i'||LP((l,oo),rdr-) + 



Furthermore, if p ^ 2, 

drW 



drW w{r) 



LP {{1,00), rdr) 



< C{p) (||/||LP((l,oo),r-dr) + |«|) 



LP {{1,00), rdr) 



+ 



LP {{1,00), rdr) 



<C{p) 



LP {{1, 00), rdr) 



+ \a\ + ep\h\) , 



(2.38) 



(2.39) 



with Ep = 1 if p > 2 and ep = if p < 2. 
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Proof. We define w = w{r)/r for r > 1. Then, w satisfies 

rdrrw{r) + 3drw{r) = f{r), forrG(l,oo); (2.40) 
drw{l) = a, (2.41) 

Following the method of the proof of Proposition 12.61 we multiply ()2.40p by \drw\^~'^ drW on [1,-R]. 
After integration by parts, this yields: 

/ f\drwF~'^drWrdr = - \r'^\drw\^]^ + ( 3 — -] [ Idrwl^rdr, 
Jl P ^ \ Pj Jl 

We conclude that: 

\\dr'w\\LP{{l,oo),rdr) < C{p) (||/||LP{(l,oo),rdr) + . (2.42) 

Expanding drW, we remark that drrW = f — drW so that ()2.42p implies ()2.38p . 
If p ^ 2, we then apply Proposition 12.71 to drW. This yields 



,,,,,,, . r ^ . ^ ^ (ll/llL.((l,oo),rdr-) + |a| ■ (2.43) 

LP((l,oo),rdr-) LP((l,oo),rdr) \\ ' n / / 

Since drW = dr-w/r — wjr^, estimates (|2.42p and (|2.43p immediately yield (|2.39p . □ 

3. Study of solutions to (ll.22l) - (ll.29p 

The ultimate goal of this section is to prove Theorem 11.11 and Theorem 11.21 In all this sec- 
tion, we assume that v = \ for simplicity. This can be done without loss of generality by setting 
(Vi/(t, x), Pi/(t, x)) := iy {t / u^x) P{t / V, x) / v) . Because of the computations we presented in the pre- 
vious section, we first analyze separately the decay of solutions to the Stokes equation with a fixed 
obstacle and then, we compute the long-time behavior of solutions to both heat equations with dy- 
namic boundary conditions. We conclude by combining all these computations. 

3.1. Decay of solutions to (12723]) - (l2:25]> . System ([223]) -([225]) has aheady been studied in the 
frame of ^^{Fq) spaces [5j Theorem 1.2], [6]: 

Theorem 3.1. For each q G (l,oo), the Stokes operator of the linear problem ()2.23p - ()2.25p generates 
a semigroup S'ij(t) on L%{Tq). Moreover, this semigroup satisfies the following decay estimates for 

VR{t,-)= SR{t)VR{0,■)■■ 
• For p £ [q,oo], there exists = Ki,R{p,q) > such that for every vr{0,-) G La{J^o), 

\\vRit,-)\\L^^iTo)<I<i,Ri~'"~'\\MO,-)\\L%{To)^ for all t > 0. (3.1) 

• If q < 2, for p G [q, 2], there exists i^2,R = K2,r{p, q) > such that for every vr{0, •) G La{To), 

\\VvRit,-)\\LP^^^^)<K2,Rt--^^-v--^\\vRiO,-)\\Ll(To)^ for all t > 0. (3.2) 

• For p G [max{2,q},oo), there exists K^^r = K^^R^p^q) > such that for every vr{0,-) G La{To), 

iiv»H(v)ii„,.,o<|^''«;l;\;||"f f"" (3.3) 

[ K3,B t " W'VBiO, ■)\\lI{To) . for all t > 1. 

For localized initial data it is possible to obtain a much sharper description of the long-time behavior 
of vr by following the spirit of our spherical-harmonic decomposition. To this end, we need a general 
result on the decay of solutions to heat equations with dynamic boundary conditions. This result is 
detailed in the following subsection. So, we postpone the more precise computation of the long-time 
behavior of vr^ to the end of this section (see Theorem 13. 13p . 
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3.2. Semigroup estimates. We proceed with the computation of the long-time behavior of solutions 
to (|2.1'2|) - (|2.14|) and (|2.28|) - (j2.3U|) . We note that both equations are examples of the family of systems: 

dty + { —drirdrv) + ) = , for (t, r) G (0, oo) x (1, r) ; (3.4) 

y{t, 1) = £y (t) , for (t, r) G (0, oo) x (1, r) ; (3.5) 
/y(t) = d(9^y(t,l) -%(*,!)), fortG(0,oo); (3.6) 

with parameters a > and /c G N U {0}. Indeed {z,iz) solution of ()2.28p - ()2.30p is a solution to 
^i^-^M in the case 



k = a 



■K + m 

whereas {w,lw) solution of ()2.12p - ()2.14p is a solution to ()3.4p ~ (|3.6p in the case 

27r 

k = 1 ^~~J' 

To compute the decay of solutions to (|3.4p - (|3.6p . we use classical methods for parabolic equations 
(see [3 [231 [TO]). In our context, due to the presence of the solid, we shall refer extensively to the 
works |19[ [IB] of A. Munnier and E. Zuazua which study thoroughly the equation 

[ dtv - A^nv = 0, for (t, x) G (0, oo) x R"\5(0, 1), 

v(t, x) = i^{t), for {t, x) G (0, oo) x 5"-^, (3 7) 



Kit) = a 



/ drv{t, x) da, for t G (0, 00) 



where a > is a fixed real number. Formally, for arbitrary fc G N, (y, iy) is a solution to ()3.4p - p.6p 
if and only if the pair (v, ly) defined by 

£y{t)=£Y{t), v(t,r,w) := Vr>l, Vw G ^'^'S (3.8) 



is a solution of equation (j3.7p for 



r 



n = 2fc + 2 , a = . (3.9) 



In this subsection, we fix /c G N U {0} and a > and study the long-time behavior of the solution 
of system (|3.4p ~ (|3.6p . By (j3.9p . this fixes also values for n and a. 



In order to study system ()3.7p . A. Munnier and E. Zuazua introduce the functional spaces 
i2P(M") = {y G LP(]R"), Vy = in 5(0,1)}, {p G [l,oo]), 
endowed with the norm: 

ll^ll£P(R«) = llyllLP{R"\B(0,l)) + ^l^^l^' "^^^^ P < 00 , 

ll'5^ll£°°(R") = max(||7y||ioo(Kn\B(o,i)), |^y|), corresponding to p = 00 , 

where iy is the mean value of Y in the ball: 

As before, in what follows, we identify (v,£v) e LP{W \ B{0, 1)) x R with the extension V G £P(M") 
given by V = 1_b(o,i)^v + 1r"\s(o,i)V, and we shall write V = (v, ^v) to denote this extension. 
We also introduce a radial variant of -£P(M^)-spaces: 

QP ■= {Y = {y,ly) radial function, such that Y G ^^'(M^)} . 
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This space is endowed with the norm: 

2-K 

||y||f;oc. = max(||y||^oo(jrg), |-^y|), corresponding to p = oo . 

In the case p = 2, this space is a Hilbert space associated with the scalar product: 

r 2'K 

For p ^ 2, extending this scalar product by a density argument enables to identify the dual of £^ with 
£P where p' is the conjugate exponent of p. 

With these notations, A. Munnier and E. Zuazua prove in [W\ I18|: 

Theorem 3.2 (Decay estimates for ([SZD, [ISllIS]). Given (vo,4o) e there exists a unique 

solution 

(v,4) G C([0,oo);£2(M")) of ([321) such that {v{0, ■),l^{0)) = (vo,4o)- This solution satisfies: 

||(v(t,-),4(t))||£2(I,n) < ||(vo,4o)k2(M"), Vt>0. (3.10) 

Moreover, i/(vo,£vo) ^ for some q G [l,oo], for all p £ [q,oo], there exists a constant C{p,q) 

such that 

tt(V.-i/p)||(v(t,.),4(t))||^,(^„)< ||(vo,4J|U,(r), Vt>l. (3.11) 

Theorem 3.3 (First term in the asymptotic expansion of solutions of p.7p . [ISllIB]). Given (vo,^vo) G 
£2(M") such that vq G ^^(M" \ S(0, 1); exp(|x| V4)dx), setting 



Jr"-\B(Q,1) O! 



-\B(0,1) 

we get 

• for allt>0 andpG [l,oo], (v(t, •), 4(i)) G £^(M") 

• /or all p G [1, oo], there exists a constant Cp such that for all t > 0, 

tt(i-iM .) _ MG(t)||^,(i,„\5(o,i)) < CpRi,p{t), 

M 



t2 

where 



{ATity. 



< CR2{t), 



1 / bp 



G{t,x) = exp 

(47rt) 2 V 4t 

ancf, denoting by 5n,2 the Kronecker symbol: 

/ (5„,2|log(t)| + l)t-i/2 z/pG[l,2], n (p-l)(p-n) 

-Ri,p(i) = < ,,,, , ,,,_i/2+fl ^ o with9n,p = -- 



(5„,2|log(t)| + l)t-i/2+^"- i/p>2, "'^ 2p(2p + n(p-l))' 

i?2(t) = «2|log(t)|l/2 + l)t-l/("+2). 

We do not give a comprehensive proof of Theorems I3.2H3.3I and let the reader refer to \18\ [T9] for 
further details. Let us only recall that the proof of Theorem 13.21 is based on the remark that ()3.7p 
reduces to the abstract ODE: dtV + ^mzV = 0, where Amz is the unbounded operator with domain 

V{Amz) = {V = (v, 4) G H\W\B{0, 1)) X M with V||,|=i = 4}, (3.12) 

such that: 



^r„.(v,4)= ^ . , . ; = f . , \ ■ (3-13) 





A. Munnier and E. Zuazua show that this operator is maximal monotone which implies the existence 
of a contraction semigroup on £^(]R'^) representing the unique solution to ()3.7|) . Further classical 
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smoothing properties of this semigroup also yield that, for (vo,^vo) ^ T^i^mz), the unique solution to 
(j3.7p satisfies: 

VGCH[0,oo);£2(]R"))nC([0,oo);P(^^,)). (3.14) 

We remark that (|3.7|) is rotational invariant. Hence, considering radial data and noting that transfor- 
mation (|3.8|) is a bi-continuous one-to-one and onto mapping from ii^ to radial functions in £^(]R") 
(with n = 2k + 2), Theorem 13.21 implies: 

Theorem 3.4. Given Yq G £^ there exists a unique solution Y G C([0,oo);£^) of ()3.4p - (|3.6p such 
that Y{0, ■) = Yq. This solution satisfies: 

\\Y{t,-)\y <\\Yo\y, yt>o. 

This theorem implies again that the solution to (j3.4p - (j3.6p is given by a contraction semigroup on 
£^ denoted by Sy in what follows. The results in [THl IS] are not sufficient for our purposes. Indeed, 
we also have to compute decay rates in il^ — 2'^ spaces, similar to the ones in (|3.1ip . But, when n ^ 2 
(equivalently k ^ 0) and p ^ 2, the transformation (j3.8p is not an isometry between £''(]R") and 2^, 
so that the "change of dimension" argument does not yield the expected result. Besides, we will also 
derive estimates on the drU, y/r, drrV, dru/r, and y/r'^ in LP{To), for which no precise estimates were 
given in [191 [18], except in the case p = 2. 

In the following subsection, we adapt the arguments of [191 H] to system (j3.4p - (j3.6p to estimate 
the decay of Sy in We then explain how to derive estimates on the derivatives of solutions of 
dSaD-D in QP. 



3.2.1. — 2'' estimates on y. Inspired in \19\ [T8] , we prove the following £^ — 2'^ decay estimates 
for solutions of (SM-<SM- 



Theorem 3.5. For all q G [l,oo), system (|3.4p - (|3.6p is well-posed in 2'^: given Yq 2'^ there is one 
unique solution Y of (|3.4p - (|3.6p in C([0, oo); £''). This solution satisfies: 

r(i)IU, < \\Yoh.. (3.15) 

We furthermore have the following 2^ estimates: for all p G [g,oo], Y belongs to C((0,oo);£^) and 
there exists a constant C such that 

ty<i-yp\\Y{t)\\^, <C\\Yo\\^,, t>0. (3.16) 

Furthermore, ifYQ also belongs to 2°°, we also have \\Y{t)\\^ao < ll^ll^oo- 

Before going into the proof of Theorem 13.51 let us emphasize that estimates (|3.15p - (|3.16p are 
different from the ones in (j3.1ip when k = 1, i.e. n = 4, that corresponds to {w,(l.w) solutions of 
(|2.12p - (j2.14p . To be more precise, in that case, using the transformation (j3.8p for r > 1, (|3.1ip would 
then read: for all q G [l,oo], p G [(?,oo], there exists a constant C such that for all Wq = {ujo,£wq) 
satisfying wo/r G Li{R'^ \ B{0, 1)), the solution W = {w,iw) of (|2T2D - (l2Ti]) satisfies, for ah t > 0, 

w{t) 



t 



2(l/g-l/p) 



< c 

Lp(M4\B{0,1)) 



V r ' 



£9 



Hence, the solution W of (j2.12p - (|2.14p will simultaneously satisfy the decay estimates (j3.16p and 
(|3.17p . Actually, as we explain below, both results can be proved following the same strategy based 
on suitable multipliers, the only difference being Sobolev's embeddings. 

Proof Let Yq G 2^ and Y = {yjy) G C{[0, oo); 2"^) be the unique solution to i^i^-i^M) given by 
Theorem 13.41 Up to assume that Yq is sufficiently smooth and vanish sufficiently rapidly at infinity 
we can apply the regularizing effect of the semigroup in M" (see ()3.14p ) so that, going back in we 
have Y G C([0, oo); i?^ (M^)) and y G C{[0,ooy, H^iJ'o)) ■ Then, the idea is to multiply equation ([TID 
by j'{y) for smooth non-decreasing convex functional j = j{y) with at most linear growth at infinity. 
After integration by parts, this yields 

i { I M + '^-^ii^y)] + / j"(2/)|VyP + ^/(^y)^y + ! ^ = 0. (3.18) 
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After a classical regularization argument, one can show that such estimate can be extended to the 
convex functionals j{y) = \y\'^, for q £ [l,oo), and this yields: 



dt 



Y\\l,) < 0. (3.19) 



Similarly, using functionals of the form j{y) = (y — K)^, after a suitable regularization argument, one 
derives 

^(imk~)<o. 

Based on the contraction property (j3.19p . the semigroup Sy{t) can be uniquely extended by density 
to initial data in as an operator from to itself. We thus have the well-posedness of (|3.4P " (|3.6p 
m any £,1, q £ [1, oo) (£2 n H^{R^ \ B{0, 1)) is not densely embedded in thus our argument does 
not apply in that case). This yields also the decay estimates (j3.15p for all q £ [l,oo] and yo £ 
Note that the decay estimates (|3.15p also coincide with the decay estimates ()3.16p for any p = q < oo. 

Actually, one can go even further. Taking j(y) = for p > 2, estimate (j3.18p implies (forgetting 
the two last terms which are non- negative) : 

±nn.uP ^ , 4(P-1) 
dt P JTo 

Using then suitable Sobolev embeddings and interpolation estimate (actually, this is the only step 
where the dimension plays a role), one gets Lemma 2.2 in [19j (the proof is done in [IH]), and in 
particular \19\ (2.17)]: there exists a constant C such that for all functions Y = {y^ly) £ with 
y G H\Fo): 

\\Y\\%.<C\\Y\\l4Vy\\l,^:,^y 
Applying it to we get the existence of a constant C such that for all g > 1, 

Plugging this estimate in (|3.20p for p = 2q and using the fact that the H'^'-norm of y decays according 

to ^7m . 



y\\%) + ^^^^l |V(|yr/2)|2<o, (3.20) 



Of course, this implies that there exists a constant C independent oi q £ [1, oo) such that 



s(ii^(')iig.) + z^(ii''<"ii?-)'^»' 



This yields the following decay property: there exists a constant C > independent of g > such 
that for all q £ [l,oo), 

t\\Y{t)\\%,<C{m\%f. (3.21) 
Then, the iteration argument of [26] based on (|3.2ip applies and yields 

t^/'^\\Y{t)\\^^ <C\\Y4^,, t>0. 

Other estimates in ()3.16p are deduced for arbitrary p £ [q, oo) by interpolating the cases p = q and 
p = CO. □ 

As we mentioned in the above proof, the semigroup Sy associated with system (|3.4p - (|3.6p extends 
to a semigroup on 2fl for all q £ [l,oo) that we still denote the same for simplicity. Consequently, 
Corollary Corollary 2.5, p. 5] implies that it is associated to a closed linear operator. In this case 
the operator reads Aq where 

V{Aq) = {Y = {yjy) G with AqY £ 
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and 

U2 



AqY = Aq{yjY) 




( . ^ y 



a 



\ I dryda + kdiv 



3.2.2. £,P — 2,'^ estimates on dtY. In the case p = 2, as A2 is self-adjoint (see [181 App. A]), Theorem 
7.7 in [2] states that, if Yq £ £^ the solution Y of <^i^-i^M) belongs to C°°((0, 00); n£eN^(^2)) and 

\\dtYit)\\^, = \\A2Yit)\\^,<j\\Yo\\^,. 

Extending this result to the case, for q € (l,oo) turns out to be slightly more intricate. 

Theorem 3.6. For all q G (l,oo), there exists a constant C = C{q) such that for all Yq £ 2fl , the 
solution y of ()3.4p - ()3.6p satisfies, for all t > 0, 

\\dtY{t)\\^,<j\\Yo\\^,. (3.22) 

Proof. The proof of such result is rather classical, but we did not find precise reference in our precise 
setting. We follow the proof of Theorem 3.6 in |20^ Chapter 7]. First, we recall that is a Banach 
space whose dual is identified with , for q' = q/{q — 1), when taking the duality pairing 



r 2tt 

(Yi,y2)£9 = / yim + —^Yi^Y2, 



for Yi = (yi,£y^), Y2 = (2/2 ) ^¥■2)- Note that, in this proof only, we extend S*' to functions having 
complex values. We focus on the case q > 2. 

For Y £ V{Ag), Y* = \ Y\''-W belongs to £«' and satisfies 

{Y,Y ) (^q j^q' — and W^i' ~ ll^ll^? • 

Besides, easy computations yield 

{A,Y,Y*),,^,=lf \y\i-^\Vy\^+fl-l) [ \yr\yVyf + k^[ ^ + 27rA:|^y j". 

In particular, both first terms can be expressed easily in terms of \y\'^~ yVy. So, we introduce the 
vectors a = a{x), and 6 = h{x) of defined by |y| 2~^yVy = a + ib. We get 

{AqY,Y*)^,^,, ={q-l) [ |ap+ / \bf + {q - 2)i [ a-b+k^[ + 2Trk\g\'' . 

In particular, 

> (g-l)||a||i.(^^,) + ||6||i2(^^), 

whereas 
This implies 



i{AqY,Y*)^,^g']\ l\q-2\ , , 

^ ^ < - , . (3.23) 



^({AgY,Y*)^,^A 2^0^ 



From Theorem 13.51 —Aq generates a Cq semigroup of contractions on hence Theorem 3.1 in [201 
Chapter 1] implies that for all A > 0, A is in the resolvent set of —Aq. 

For g > 2, from (|3.23|) . the numerical range S{—Aq) is contained in the sector Egg = {AgC\{0} : 
|argA| > vr — 6q} where 

00 = arctan G [0,vr/2). 
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In particular, choosing 6i G {6q,tt/2), denoting by = {A G C \ {0} : |argA| > vr — ^1} the 
corresponding sector of C, and using the fact that is in the resolvent set, Theorem 3.9 in [2^ 
Chapter 1] implies the existence of a Cg = CQ{q) such that 

||(A/ + ^,)-i||^(.p)<^, VAGC\Se,. (3.24) 

Now, the regularizing properties of the semigroup generated by —Aq are a consequence of Theorem 
5.2 in j20| Chapter 2] and the above resolvent estimate. However, here again, we need to be careful 
since Theorem 5.2 in Chapter 2] requires that belongs to the resolvent set of Aq, which is not 
the case here. Set 62 G (7r/2, vr — 9i). For each e > 0, we introduce the curve F^, defined for e > by 
the path composed as follows: 



F, 



-yoexp(-i6'2), 

eexp(i0), 

/9exp(i6'2), 



6 G 
P e (e,ooj. 



00, —e). 
^2, O2), 



oriented in the increasing directions of the parameters. Then, for t > 0, we use the formula 

Syit) = ^J^ e'\XI + Aqr'dX. 

This integral converges due to the resolvent estimates (j3.24p and can be differentiated with respect to 
time since 

1 



2tt 



-5R(A)t 



X{\I + Aq)-'\\^^^^^^ dX<^ + Cte, 



where the constant Cg^ does not depend on t and e > and the constant Ct depends on t but not on 
e > 0. Of course, letting then e — )• 0, this yields 

C92 



dX 



\dtSy 



< 



t 



This completes the proof of (|3.22p for g > 2. The case q G (1,2) can be deduced by a simple duality 
argument. □ 

Remark 3.7. Actually, following the proof of Theorem 5.2 in [20\ Chapter 2], one can prove that —Aq 
generates an analytic semigroup on £^ for all q G (1, 00). 

In the two next subsections, we apply the semigroup estimates we have proved to systems ()2.28p - 
((OnD and (l2l^ - (l2TiD . 

3.3. Decay of solutions to (127281) (l2T30l> . We first consider the solution Z = {z,iz) of ([2:28]) - 
(|2.30p . As we noticed previously, this corresponds to the computations of the previous subsection in 
the case 

47r 

k = , a = . 

TT + m 

We obtain in this way the following decay estimates on solutions: 

Theorem 3.8. Given q G (l,oo) and radial Zq G -G*^, there exists a unique solution Z G C([0, 00); H"^) 
to (|2.28p - ()2.30p such that Z(0, •) = Zq. This solution satisfies the further decay estimates: 

• for all p G [q,oo] we have Z G C((0,oo);£^) and there exists a constant Ki^i = Ki^i(p,q) such 
that: 

\\Z{tr)Uv <Ki^itl-\\\Z4st'^, Vt>0, (3.25) 

• if q < 2 for all p G [q,2), we have {drZ,z/r) G C{{0,oo); IJ'{To)) and there exists i^2,i = 
K2,i{p,q) such that: 



\drZ{t, •) 



+ 



< K2 I t ^ P 1 ||Zo||£<J , 



LP (To) 



(3.26) 
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if q £ (1, oo), for all p G [max{2, q},oo) with p > 2 we have {drZ,z/r) G C((0, oo); L^(-Fo)) o,nd 
there exists K^^i = Ks^i{p,q) such that: 



\drZ{t,-) 





z{t,-) 




+ 








r 





Ks^it 2 + P ^||Zo||£9 , Vt < 1, 
1 1 ||Zo||£9 ) Vi > 1 . 



(3.27) 



These decay estimates are also satisfied for g = 1 and p G (1, 00) \ {2}. 



Proof. Existence of solutions and £P — Sfl decay estimates are straightforward applications of Theorem 
13.51 in the case k = Q and a = 47r/ {-n + m) > 0. We focus on estimates (|3.26p - (|3.27p . Actually, we only 
need to prove the case p = q, as other cases are then obtained by combining the estimates (j3.26p - 
(j3.27p for p = q ^ 2 between t/2 and t with (j3.25p between and t/2. Indeed it will follow from the 
semigroup property: 



\drZ{t,-)\\ 



+ 



z{t, 



LP {To) 



<K2,i(j},p) {t/2)-2\\Zit/2,-)hp 

< K2,i{p,p)Ki^i{p,q) {t/2y'^+p-'<\\ZoU,. 



For radial G -C, estimate ()3.22p implies 

drZ{t) 



drrZ{t) 



LI (To) 



(J 

+ \drZ{t,l)\ < -\\Zo 



We are now in position to apply Proposition 12.61 to drZ which yields that (see ()2.36p ) 

drz{t) 



'rrZ{b)\\Li{To) 



+ 



<^IIZII 
S y||^o||£9 



(3.28) 



From the Gagliardo Niremberg inequality in exterior domains, see [1], we have then: for q G [l,oo], 
for all z such that dxxZ, dxyZ and dyyZ belong to L'^{To), 



(3.29) 



Since we are focusing on the case of radial solutions, estimates ()3.28p - (|3.29p and the fact that for 
radial functions 



\\dxxZ\\Li{To) + \\9xyZ\\Li{To) + ll'^TO^IIi'' W)) - ^ W^rr z\\ Li (Tq) 



dr-Z 




r 


LiiTo)) 



imply 



C 

\drZit)\\L',(J^,,) < -^||Zo||£9. 



To conclude the proof of Theorem 13.81 we prove the boundedness of the mapping drZ 1— )■ z/r. As 
z G L^{J-q), this is already contained in Proposition 12. 7| provided we get a suitable estimate on 
z{t, 1) = izit)- But, using (|3.25p for p = 00, we get 



Thus, (IOTP implies: 



\izit)\ = \zit,l)\ <Cgt—i\\ZoUi. 
z{t) 



Li (To) 



0||£'J 



where = 1 if g > 2 and = if g < 2. We obtain ()3.26p and (|3.27p comparing the size of the 
different terms on the right-hand side depending on g < 2 or g > 2 and t>lort<l. □ 
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3.4. Decay of solutions to (l2n^ - (f27T4l> . The equation (f2J2]l - (f23i]l of w is linked to the com- 
putations in Section in the case k = 1 and a = J'/2it. Thus, we compute the following time-decay 
of solutions: 

Theorem 3.9. Given q G (1, oo) and radial Wq G il^, there exists a unique solution W G C([0, oo); 
to (|2.12p - (j2.14p such that W{0, •) = Wq. This solution satisfies the further decay estimates: 

• for all p £ [q,oo] we have W G C((0,oo);£p) and there exists a constant Ki^ = Kifi{p,q) such 
that: 

\\W{t,-)hp <Kifit-p--^\\WoU,, Vt>0, (3.30) 

• if q < 2 for all p G [(?, 2), we have {drW,w/r) G C((0, oo); LP(Jx))) o,nd there exists ^2,0 = 
K2fi{p.,q) such that: 

W{t,' 



\drW{t,-)\\LPiJ^o) + 



<K2,ot ^||VFo||£5, Vt>0. 



(3.31) 



ifq G (1, oo), for allp G [max{2, q}, oo) satisfying p > 2 we have {drW,w/r) G C((0, oo); L^^Fq)) 
and there exists K^^ = K3fi{p,q) such that: 



\\drW{t,-)\\Lp(^^^) + 



wit,-) 


4 


r 


LP [To) { 



K^fit—^^^v—^WWoU'^, Vt<l, 
i^3,ot"^||VFo||£, , Vt>l. 



(3.32) 



Proof. Again, existence of solutions and — £fl decay estimates are straightforward applications of 
Theorem 13.51 in the case k = 1 and a = J /2tt > 0. We focus now on gradient estimates in the case 
p = q ^2, the estimates (|3.3ip - ()3.32p with p > q, p ^ 2 being a simple consequence of the semigroup 
property. 

Let q G (l,oo). We note that w G Ci((0, oo); L'?(J"o)) with \\dtW{t)\\s,q < C/t\\WoUi yields 



drrwit) + - ^ G Cmoo);LHTo)) 



drW{t) w{t) 



drrW{t) + 

"LI {To) 

Recalling estimates and (ITMD . for all t > 0, 



C. 



+ \drWit,l)-wit,l)\ < y||l^o||£^ 



\\drrW{t)\\L,iT,) < y||M^0||£., 



drW{t) 



Li{{l,oo),rdr) 



+ 



wit) 



^ ^ <-\\Woh^+Ciq)eg\wit^)l 

Li{{l,oo),rdr) t 



with £q = 1 ii q > 2 and Eg = if q < 2. But, for g > 2, estimate (|3.17p with p = oo yields 



\£wit)\ = \wit,l)\<Ct-^/l\\(^^,iWo) 
where the last estimate is a consequence of g > 2. Hence 



< C7t~2/9||Wf 



0||£9, 



dr-wit) 



L'J{(l,oo),rdr) 



+ 



wit) 



L9((l,oo),rdr) 

We can then bound ||(9a;a;tt;||/,q(jr^), ||i9a;j/t(;||ig(jrjj) and ||(9yyi(;||ig(j-Q) in the same way as in the previous 
proof. Applying interpolation inequality ()3.29p to W we then obtain that drW G C((0, oo); L'?(/"o)) 
with: 

||a,u;(t)||i,(^„) < C {t-^i^ + Egt-^i'^) \\W4.,. 
To get the decay of w/r, we then simply use that 

wit) 2 wit) 



Li{{l,oo),rdr) 



< 



Li{{l,oo),rdr) 



\\wit)\\ 



21- 



□ 
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3.5. Decay estimates of solutions to the Stokes system. It remains now to combine together 
the results obtained in Subsections I3.H 13.31 and 13.41 to prove our main results regarding the long-time 
behavior of Stokes solutions. 

3.5.1. Proof of Theorem 

Proof. Given q £ (l,oo), as 

C^°°(M2) n C'^ is dense in we remark that it is sufficient to prove decay 
estimate for initial data Vq £ C'^ Ci C^(M^). We emphasize Vq S for all q G [l,oo] under these 
assumptions. We denote (Wq, ^q,'^o, Vr^q) the spherical harmonic decomposition of this initial data. 
We already know that there exists a unique solution to (fr22]) - (fr29]l in C([0, (x>);C^) for such an initial 
data. 



First, we decompose the solution S{t)VQ of (|1.22|) - (jl.29p into the spherical-harmonic decomposition 
{W, ^, Vr) of Proposition 12. 1[ According to Proposition 12.21 this decomposition satisfies: 

• We C([0, oo); £2) and is a solution to (f2l^ - (im] ): 

• VrG C{[0,oo);Ll{To)) and is a solution to ([223]) -([225]). 

According to Theorem 13.91 and Theorem 13. H these are the respective unique solutions to ()2.12p - ()2.14p 
and ()2.23p - ()2.25p in these spaces, with respective initial data Wq and vr^. We can then apply the 
decay estimates of Theorem 13.91 and Theorem 13.11 to these solutions. 

Referring moreover to Proposition 12.31 and Remark 12.41 we have: 

• Z$ E C([0, oo); £2) and is a solution to ([228]) -([230]); 

• £ C([0,oo);£2) and is a solution to ([228]) -([220]). 

Consequently, applying Theorem 13.81 these are the unique solutions to ()2.28p - ()2.30p in these spaces, 
with respective initial data Z^^ = —{dj-^Q + ^o/r) and Z^^ = dr'^o + '^o/r and the decay estimates 
of Theorem 13.81 are also satisfied by Z<j, and Z^ . 

We proceed with — L'^ estimates. Applying ()2.3p we have: 

l|W^o||£9 + ll^-i'.olU'j + ||-2'*,o||£9 + II^J?,o|lL«(Jb) - ^gll^ollo • 
Then, combining decay estimates of the different components in the spherical-harmonic decomposition 
obtained in ([21]), (IX^ . and (^150]) . we have for t > 0: 



On the left hand side, we have for instance Z^{t) = dr'^{t) + ^{t)/r with ^'(t, 1) = Z^^{t^ l)/2 so that: 



|^'(t,i)| <c\\z<,{t)Uv 



Hence, applying Proposition 12.61 we obtain: 

l|f^r^||LP{(0,oo);rdr) 



+ 



LP{(0,oo);rdr) 



<C||Z^||£P. 



Applying similar argument to bound we finally obtain that: 



||^^(t)IUp + ||V^fiWllLS(.Fo) + 



LP((0,oo);rdr) 



+ 



LP((0,oo);rdr) 



< CqtT 



CI 



\dr^{t)\\Lv{{0,oo);rdr) + 

LP({0,oo);rdr) 

Noting that ||VF(t)||£P is equivalent to ||VF(t)||LP((o,oo),rdr)) we apply (|2.4p and conclude immediately 



\\V{t)\W<C, 



1 

tp' 



CI 



We now proceed with the gradient estimates. Let us recall that the case p 
consequence of [52] and the previous inequality: 



2 is a straightforward 



\S/v 



,)<Ct-'/^Vit/2)\\c2<c/^--^-"^ 



Vr 



0\\Ci 



So we focus on the case q E (1, oo) and p E [q, oo) with p / 2. 
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Similarly to the previous computations, the method is then an application of Proposition 12.11 and 
the decay estimates obtained in Theorems 13.11 13.81 and 13. 9[ The only difference we detail now is the 
computation of 



I a 



rrV\\LP{{l,Qo),rdr) 



+ 



+ 



LP((l,oo);rdr) 



LP({l,oo);rdr) 

from the estimates for 115^2:^ ||i;,p(_7rg) and ||-2i/'/^lli''(-^o) given in (|3.'26p and (|3.27|) . We focus on ifj, 
the problem being completely similar for ip. Differentiating the definition of z^, we remark that ip 
satisfies: 

dri> _ Pi 

0. 



drrtp + 



r 

- ^(1) 



for r G (1, oo) 



Consequently, we apply Proposition 12.81 which yields: 



|9rrV'llLP((l,cx)),rdr) + 



On the other hand, we have, by definition of z^, 



LP((l,oo),rdr) 



< C\\drZ^\\ip(^jr^^-) . 



drip 



< C 



LP((l,oo),r-dr) 



r 



LP [To) 



So that, we finally get: 



15, 



'rrW\\LP({l,oo),rdr) 



+ 



drip 



+ 



LP((l,oo);rdr) 



LP((l,oo);rdr) 



< C \\drZ^\\LP(^jr^) + 







r 





This ends the proof of Theorem 11.1 



□ 



3.5.2. Duality decay estimates. For later use, based on Theorem 11.11 we derive here additional esti- 
mates on the behavior of the semigroup corresponding to (|1.22p - (|1.29|) : 



Corollary 3.10. Assume that 1 < q < p < oo and let F G Li{R'^; M2x2 
The following decay estimates hold true: 

• ^/2<(7<p<oo, there exists K4, = K4^{p, q) > such that: 



satisfying F = on Bq. 



||5(t)Pdiv F\\cP < K4{i^t)~^^p~^\\F\\Lc 



for all t > 0. 

if ^ < q <p and q <2, there exists = K^{p, q) > such that: 

K5{i^tr^^+p--^\\F\\L. 



||5(t)Pdiv F\\cP < < 



K5iiyt)-^^p\\F\\L, 



for all <t < -, 

V 

for all t > —. 

V 



(3.33) 



(3.34) 



In this corollary the divergence div is computed along rows of the matrix F . 

Before going into the proof, let us emphasize that, in our case VS{t) is not the dual operator of 
S'(t)Pdiv . Indeed if F is smooth with compact support, there holds, for all G £^ H C^{M?): 

{VS{t)ip,F) := — [ VS{t)^:F+[ VS{t)<l> : F 

\-—\ [ S{t)<^-Fnda + {S{t)<l>,dwF) 
^ ^ JdBo 

1 / S{t)<^-F7ida + {^,S{t)FdivF). 

^ ' JdBo 

Hence Corollarv 13.101 onlv concerns the restriction of the dual of VS'(t) to functions F which vanish 
at the boundary. 
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Proof. The following proof contains a construction of the operator S'(t)Pdiv on the closed subset of 
L'^(M?; M2x2{^)) of functions vanishing on Bq. We prove our result in the case 2 < q < p < oo only. 
The other cases can be done similarly. 

Let F G L'i{M?; M2x2(I^^)) such that F = on Bq. Up to a regularizing argument, we assume that 
F e C^(J"o;Af2x2(K))- Then, V{t) := 5(t)PdivF e C((0,oo);£p) for all p G (l,oo) by a straightfor- 
ward application of Theorem 11.11 For all t > and V ^ C? r\ C^(M^), we have, as S is self-adjoint 
with respect to the scalar product (•, •) we introduced on C? (see ()1.32p ): 



^,,,,< = (PdivF,5(t)y), 
divF-5(i)F, 



F ■ VS{t)V 



(as F vanishes on i?o) > 
(as Fn vanishes on dB^) 



Finally, we obtain: 



CV^CP' - ll^llL''(R2)||V5(t)y||^,/(jp.^) , 

where we apply decay estimates we obtained in Theorem ll.il as p' < q' < 2 we have from ()1.35p 



t 2 + P q\\V\ 



CP 



So that, we obtain: 

As C^' n C^(M?) is dense in C^' , this inequality implies by duality that V{t) G with norm lower 
than Ct~5+^-| ||-^||l'/(r2) . □ 

In the previous corollary we restrict p to finite values. In the case p = oo, we do not obtain a control 
of the whole solution. Nevertheless, we can obtain a result that would correspond to the case p = oo 
in ()3.33p for the translation speed iv{t) ■ This result is a new application of the added mass effect and 
relies on the fact that Kirchoff potentials are easily computed in our case. 



Corollary 3.11. Let q G [2, oo) and F G L«(]R2;M2x2 
decay estimate holds true for V{t) := S'(i)PdivF; 



\ivit)\ < Ke{q){uty 
where K^{q) depends only on q. 



2^ q 



\F\ 



LI 



satisfying F = on Bq, The following 
yt>0 (3.35) 



Proof. Let the assumptions of the corollary be satisfied. At first, we recall that we have V{t) G 
for all t G (0, oo) as has been shown in the previous corollary. We show how to prove that the first 
component iv,i of ^y(t) satisfies (|3.35p . Similar estimate for the other component iv,2 is obtained 
applying comparable arguments. 

Let ip G C°°{J^o) be given in polar coordinates by: 

cos(^) 



ij{r,e) 



y{r,9) G (l,oo) X (-7r,7r). 



Given t > we note that V :- 
of J^o- This yields: 



V{t) = {v{t),£vi't)) is divergence free on any subdomain B{0,R) \ Bq 



Idi 



V ■ VV' 



B{0,B)\Bo 



V ■ ntl^da + / V ■ UTpda 
ldB{0,R) JdBo 

Letting R — t- oo, we obtain (the exterior boundary term vanishes as f G C): 

I-2TT 



Vr(l,0)cos(0)de 



V ■ Vip ■ 



(3.36) 



28 S. ERVEDOZA, M. HILLAIRET & C. LACAVE 

We observe then that, on the one hand, we have V-i/J = \/'^(p where 

^i^r,e) = ^-^, V(r,0)G(l,oo)x(-^,7r), 
r 

on the other hand: 

Vr{t,l,9)cos{0)de = 7riv,i{t)- 

JO 

Setting finally: 

S := IjtqVV^ - Ifioei 
we have that E G for arbitrary p > 1 and that (|3.36p reads: 



-{7r + m)lv,i{t)=miv{t)-iE+ v{t) ■ ^ = {V{t),E) , Vt>0. 

Given this identity, we reproduce the computations done in the proof of the previous corollary. We 
obtain: 

- (vr + m)iv,iit) = I F - yS{t)E , (3.37) 

which implies 

The proof now reduces to find a bound on ||VS'(t)H||j^,/j-j-^^. To this end, we remark that the 
spherical-harmonic decomposition of H reduces to the first mode (p{r,9) = min{r, 1/r} sin(^). Going 
back to the computation of Section 2, we note that S{t)E is given by its first-mode, corresponding to 
= —dr^ — ^/r where $(r) = min{r, 1/r}. This mode satisfies (|2.28p - ()2.30p with initial condition: 

■^<i>,o = —2 • 1b„ ■ 

Consequently, for q > 2, Z^^ G and we apply Theorem 13.81 with "p"= q' and "g"= 1, which yields 
(see (1X26]) ') 



.Z$(t,-)|lL,'(j-o) + 



Z$(t, ■ 



<K{vt) 2 + ^ = (i/t)"2-^ ^ Vt>0. 

Li' {To) 



r 

We go back to E as in the proof of Theorem 11.11 so that: 

\\VS{t)E\\^,,^^^^<K{vt)-'--^>, Vt>0. 

Plugging this estimate in ()3.37p we obtain the expected result for q > 2. 

The case q = 2, corresponding to q' = 2, does not immediately follows from Theorem l3.81 but rather 
from the fact that 

||V^ fcD (l,r] sm{e)] \y < C\\Z^ f^,-^ 11^2 < C{vt)~^'^\\Z^fl\\^. 



and from the £^ decay estimates on the gradient obtained in Theorem ll.lt 

\\VS{t/2) (v^ (c^, Q,,^ sin(^))) ||^2(^„) < C{vtr''^\\V^ Q,r) sin^) 



□ 



3.6. Asymptotic expansion of solutions to the Stokes system. This section aims at proving 
Theorem 11.21 We first show that the solutions W and Vr corresponding to the modes k ^ 1 decay 
faster than the modes corresponding to /c = 1. In a second step, we derive precisely the first-order in 
the long-time behavior of this first mode. 
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3.6.1. Faster decay on W . 

Theorem 3.12. Given a radial Wq G £^nL'^(U?,exp{\x\'^/A)dx), the unique solution to ([232]) - (fZTIl) 
satisfies W{t) S £^ for allt > and p £ [l,oo]. Furthermore, for all p G [2,oo], there exists a constant 
Cp > such that 

\\W{t)Up <Cpt-p--^\\Wor\y, 
and there exists a constant C such that 

\iwit)\ <Ct-^Wor\\i^i. 

Proof. Let us first remark that Wq G -2^ n L^(M^, exp(|j;|^/4)dx) obviously implies that Wor e 

We first focus on the decay estimate, as W{t) G £^ for i > is obvious. Given t > and p > 2 we 
apply ()3.30p with q = 2 between t/2 and t and then, we apply ()3.17|) between and t/2: 



\\W{t)U, < C7ti-^||H^(t/2)||£2 < Ctv^^ 



Wo 



<Ctp ^ ||VI7Qr||£i . 



Concerning the second estimate, it suffices to use (|1.3ip and ()3.17p : 



W{t) 



< cr 



Wo 



£1 



= Ct-^\\Wor\\i^i. 
(see Proposition 12. ip verifies (jl.SSp . 



□ 



\^w(t)\ < 

This ends the proof. 

As a consequence, we can already note that u}s(t)Vo 
3.6.2. Faster decay on Vr. 

Theorem 3.13. Given Vufi G Ll{J^o) n L^(]R^, exp(|xp/4)da;); for all p G [2,oo], there exists a 
constant C = C{p,viifi) such that: 

Proof. In order to prove Theorem 13. 131 we expand Vr solution of (|2.23p - ()2.25p on its Fourier basis: 



VRit,r,e) = V' 



r cos 



+ Lpk{t,r)sm.{ke)) 



k>2 



where ipk{t, 1) = drip^it^ 1) = Vkit^ 1) = drfkit, 1) = thanks to the homogeneous Dirichlet boundary 
conditions satisfied by the restriction vr of Vr on J^o- Note that vr does not contain any or 1 mode 
due to the orthogonality condition (|2.2p . 

As in the case k = 0,1, we can show that for all A; > 2 the new unknown = ■ = 

l/r^dr[r^'4)k{t,r)\ or Zk = z^^k{t,r) = -l/r''dr[r''(pk{t,r)] satisfies: 



dtZk + ( --dr{rdrZk) + ^^—^P-Zk 








for (t,r) G (0,oo) x (l,oo); 
for t G (0,oo). 



Zk{t, 1 

One can then use the asymptotic formula given by Theorem 13.31 for 

VR{t, r,9) = Y^ iz^,kit, r) cos((A; - 1)9) + r) sin((/c - 1)9)) 

k>2 

which is a solution of ()3.7p for n = 2, arbitrary q > and vanishing initial mass M = 0. This 
immediately yields that, provided 

^^(O) G L2(Jo,exp(|x|V4)dx), 
which holds true since Vr^ is assumed to belong to L^(M^, exp(|xp/4)dx), we have 

t'-'^n\Mt,-)\\LnT,^<CRi,p{t). 
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In particular, for p = 2, this implies that 



/oo 
r(|z^,,.(t,r)p + |z^,fc(t,r)|2 dr < CRi,2{tf 



But recall that, for k > 2, 



and V'fc(^) 1) = 0. Hence for all R > 1, 



R 



\z-4,,kf' rdr 



R 



\dri^k\^rdr + 



R 



As ^fc(l) = 0, passing to the limit ii —)• oo, we get 



\dri>k\'^rdr + k"^ 



rdr < 



rdr + k i (?r (IV'fcP) dr. 



and thus. 

Using then the semigroup estimates we get, for p > 2, 

This concludes the proof of Theorem I3.13t as Vr simply vanishes in i?(0, 1). 



□ 



3.6.3. Proof of Theorem ] l.S\ 



Proof. Let Vq ^ r\ L2(m2^ exp(|x| V4)dj;) and V{t) the unique associated solution to ^T?m - ^T?M . 
Note that Vq G for all q G (1, oo) so that we already know that V{t) G D' for all p G [2, oo) for all 
t > from Theorem ll.il Let now p G [2, oo) and (VF, Vr) the spherical-harmonic decomposition 

of y. 

The components ^ and W are computed as means of y in so that they inherit the asymptotic 
decay in r of the data Vq. Combining this remark with Proposition 12.11 this yields that: 

1^(0,-) G £^nL2(M2^exp(|x|V4)drc), Vr(0, •) G (Jg) n L2(m2, exp(|x| V4)drc) 

(Zci>(0,-),^*(0,-)) e -C^nL2(M2,exp(|3;|V4)drc). 



Consequently, Theorems 13.131 and 13.121 imply respectively: 

log(i)l 



\\yR{tr)\\Li(^T-^) 



O 



\\w{t,-)hv = o{t'/p--'i^) 



j3/2-l/p^ 

We focus now on Z$ and Z^jJ. Using Theorem 13.31 with a = 2/(7r + m), we immediately get: 



t^-^/P ||z$(t, •) - M$G(t)||ip(^,,) < Cpi2i,p(i), t 
t^~^'P \\z^{t, ■) - M^G(t)||^,(_^^) < CpRi,p{t), t 



Iz it) - ^ 



< CR2it), 

< CR2{t), 



(3.38) 
(3.39) 



with G and (i?i,p, -R2) as given in Theorem 13.31 in the case n = 2 and 



/oo /■! 
2'$(0,r)rdr + (vr + m) y 2'$(0,r)rdr, 

/CO /■! 
Z5,(0,r)rdr + (7r + m) y Zvt(0,r)rdr. 
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Recalling that 

Z^(t,-) = --dr(r^(t,-)) for (t,r) G [0, oo) x (0,oo); 
r 

Z^{t,r)=iz^{t) = -2$(i,r)/r = 2£i(t) for (t, r) G (0, oo) x (0,1), 
$(0, 1) = -^i(O), (by the continuity of <^>) 

and using $(0, •)/r, 9r<l*(0, •) G n L^((0, oo), exp(|rp/4)rdr), which implies the existence of a se- 
quence Rn — )■ oo such that ii„<l>(0, i?„) goes to as n — )• oo, 

M$ = (vr - 7n)$(0, 1) = (m - 7r)^i(0). 

Similarly, 

= (m - 7r)^'(0, 1) = (m - 7r)^2(0). 

By Proposition O, we recah for r G (0, 1) that ^5(4)^0,1 = -$(i,r)/r = £z^(.t)/2. From (fOHjl - ffOOD 
and the previous formulas, we have obtained (|1.37|) and (|1.39p . 

Solving $ and in terms of Z^{t, •) = {z^{t, •),iz^i't)) and Z,i,(t, •) = (zii-(t, •),£z,s,(t)), we are then 
led to define '^{t, r) on t > 0, r G (0, 00) as the extension of i{j solution of 



by 



^dr{r'4){t,r)) = G{t,r) for (t, r) G (0, 00) x (1, 00), ip{t,l) = fortG(0,oo). 



^{t, r) = — for {t, r) G (0, cx)) x (0, 1). 
ovrt 



Note that this function can be computed explicitly: 



r j_ 

27rr 
r 



1. Svrt 

Using Proposition 12. 6^ we get for all p > 1, 

i;{t,-)-M^^{t,-) 



\\dr^l^{tr) - M^drij{t,-)\\LP{To) + 
dr(p{t, •) - M^dr^{t, ■)\\lp{To) + 



-99(t,-)-M$l'(t,-) 



for (t, r) G (0,00) X (l,oo), 
for (t,r) G (0,00) X (0,1), 

< cJRi^j,{t)t^'P-^ + R2{i)t-'^), 



< Cp(Ri,p{t)t^'P-^ + R2{t)r^ 



With the expression of Ri^p and R2, we can check that —\ + ^2,p + p~l>~i~l ^r all p G [2, 00]. 
Hence for t > 1, we have: 



J'^^^^ .) ^Qg(g) _ ^^^(^^ ^^g(^) j 11^^^^^^ < 2Cpi?l,p(t), 

(^^(^^ .) ^ ^^^(^^ .) 11^^^^^^ < 2Cpi?i,p(t). 

Remark then that, denoting 



■0(t,r) 



1 



27rr 



1 — exp 



r 



< 



for (t, r) G (0,00) X (l,oo), 

we have for all p G (1, 00], 

||V^((V;(t,r) -^(i,r))cos(0))|Up(^„) < C 
We then obtain 

_ [(m - 7r)V;(t, •)(^2(0) cos(^) - £i(0) sin(^))] < CpRi,p{t). 

This yields the expected result. 



7^ 



□ 
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4. Long-time behavior of solutions to the Navier-Stokes problem 

In this section, we prove Theorem 11.31 and Theorem 11.41 We first apply Kato's method |14j of 
successive approximations yielding decay estimates for initial data Vq £ C'^. In a second subsection, 
we then extend these estimates to the case of initial data Vq G with q G (1,2] in order to get 
Theorem 11.31 We finally explain how a bootstrap argument yields Theorem II. 4i 

To simplify notations, we replace the constants Ki{p,q), Ki{p,q) and Ki{q) defined respectively in 
(fLSSD . (1333!) and (l3:35|) by Ki{p,q)v^lP~^li , K4p,q)u~^/^+^/P~^/'' , and Ki{q)u~^/^~^/i , so that the 
viscosity parameter will not appear in our computations. 

4.1. £P decay estimates for initial data. We recall that we transferred our system in the body 
frame applying the change of variable (jl.lOp . So, the equations ()l.ip - ()1.8p became (jl.lip - ()1.18p . Our 
first proposition reads: 

Proposition 4.1. Let Vq £ C"^. There exists Aq > such that, if 

\m\c^ < Ao, (4.1) 

then, the unique global weak solution V of ()l.lip - ()1.18p satisfies the following: for allp G [2,oo), there 
exists constants H{p, Aq) and //^(Ao) such that 

snpt^-p\\V{t)\\cP < H{p,Xo), and snpt^\£v{t)\ < He{Xo), 
t>o t>o 

Proof. We split the proof of Proposition 14.11 into six steps. 

Step 1: integral formulation. Following [22], we rewrite the Navier-Stokes equations (|l.lip - 
p.l8p in the following abstract form: 

dtV + AV = FF 

where 

[ev-V)-VV on To 
on Bq, 

P denotes the continuous projector from to C^, and ly is defined for V £ hy (jl.3ip . Then, 
Duhamel formula gives the following integral formulation of the above equations: 

V{t) = S{t)Vo + [ S{t- s)FF{V{s)) ds. (4.2) 
Jo 

T. Kato suggests to construct a solution by successive approximations: let the sequence (yn)eN be 
defined by 

{yI Z Smo + ICYn-u ynen, -^-^ ICYit) = I^Sit-s)FFiY)is)ds. (4.3) 

Our aim is to prove that this sequence satisfies uniformly estimates of Proposition 14.11 and converges 
for small initial data. To simplify notations, in the following we set £y„ = £n- 

Concerning the nonlinear term, we note that, ¥F{Y) is well-defined as soon as Y € for p > 2 
satisfies \7Y G LP'{Fq). Indeed, we can then split F{Y)\jr^ = —Y • VY + ly • Vy, the first term being 
in L'^(J-o) (where q = 2p/(2 + p)) and the second one in Lp'{Fo). We have then : 

PF(y) = -Pg[Vo^ • vy] + P2[i.Fo^y • vy]. 

Furthermore, we remark that, if y G (with pQ G [l,oo)) satisfies y G H^{Fq) then: 



F{V) 



F{Y) = divF(y) where F{Y) 



^Y-Y)®Y on Jo , 
on Bo . 



This property is satisfied since F(y)n vanishes on SBq as -Bo is a disk. The operator K, can then be 
defined indifferently as: 

I S{t- s)[¥F{Y){s)\ds or / [5(t - s)Pdiv]F(y)(s) ds , 
Jo Jo 
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where S{t — s)Pdiv is defined by duality. In order to get uniform estimates on tlie functions Yn and 
their hmit, we work with the second form (Step 2 to 5). In Step 6, we apply the first form to prove 
that our construction coincides with the unique global weak solution constructed in 



Step 2: estimates of t» ||y„||£8, ||y„||£2, t2 The goal of this step is to show the following 

Lemma: 

Lemma 4.2. There exists a constant Aq > such that for all Vq e satisfying ()4.ip there exists 
Ho > such that: 

3 1 
SUp{t»\\Yn{t)\\cs} < fio, SUp{||y„(t)||£2} < /io , sup 1 2 \in{t)\ < fiQ. 
t>0 t>0 t>0 

Besides, fiQ can be chosen arbitrary small, independent ofVo, up to restrict the size o/Aq- 
Proof. We are going to find by induction a sequence G„ such that for all n, 

snp{t^-^\\Yn{t)\\cs} = snptl\\Yn{t)\\cs < Gn, (4.4) 

i>0 t>0 

sup{i^"5||y„(t)||^2} = sup||y„(t)||£2 < Gn. (4.5) 
t>0 t>0 

sup{i5|£„(t)|} < Gn. (4.6) 

i>0 

It is clear from (|1.33p that (|4.4p - (|4.6|) is verified for Yq, where 

Go = max{i^i(8,2),Ki(2,2),ifi(oo,2)}||l/o||£2. (4.7) 
In the sequel, we denote by Cq the following positive constants: 

Co :=max(if4(8,4),K4(2,2),i^,(4)) (1 - r)-lT"t dr, 

where is the constant in ()3.33p and Kg in ()3.35p . 

Next, we assume that the properties are true for the rank n, and we show it for rank n + 1: using 
()3.33p with p = 8, q = 4, we get 



ts\\Yn+i{t)\\cs < Go + tsKi{8,4) / {t-s)-s\\Yn{s)\\lsds + tsKi{8,4:) / (t-s)-8|4(s)|||y„(s)||£4 ds. 

Jo Jo 

By interpolation, we have: 

ll^llr* < liynll/2'^11^11/8'^ • (4.8) 

So, we use that: 

\\Yn{s)\\ls < {s-ICnf, \\yn{s)\\c^ < S'^Gn, |^„(s)| < S'^^Gn , 

to get 

ft ft 
ti||y„H_i(t)||£8 < Go + tlKi{8,4:) / (t - s)"ts"t|Gnpds + tiK4(8,4) / {t - s^h'^Gnl'^ ds 

Jo Jo 

< Go + 2|G„|2K4(8,4) / (l-r)-8r-4dr 

Jo 

< Go + 2Co|G„p. 

Writing the same computation and using (|3.33p with p = q = 2 gives 

f-t r-t 

||y„+i(t)||£2 <Go + i^4(2,2) / {t-s)-"2\\Yn{s)\\%ds + K42,2) / (t - s)-2|4(s)|||y„(s)||£2 ds 

Jo Jo 
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Thus, we obtain 

\\Yn+i{t)\\c2 < Go + 2K^{2,2) {t-s)-h-^Gn\^ds 

Jo 

< Go + 2Co\G 

n I • 

Finally, we apply Corollary 13.111 with g = 4, which yields: 

thln+i{t)\ < Go + t^Ke{4) [\t-s)-^\Yn{s)\\lsds + t^Ke{A) [\t - s)-^in{s)\\\Yn{s)\\c^ ds 

Jo Jo 

< Go + 2Co|G'„| . 
Hence, we can take 

Gn+i = Go + 2Co|G„|^, 

in (|4.4p - (|4.6|) with Gq given by (j4.7p . Choosing Aq such that Go < l/(8Go), we easily get by induction 
that for all n G N 

O„< '-";f°^°'^ :M0. (4..) 

4C-0 

Therefore, (G„) is bounded by ^uq which implies that (|4.4|) - (j4.6p are uniform estimates. This ends the 
proof of Lemma 14.21 According to (|4.7p - (|4.9p . fiQ can be chosen arbitrarily small by taking Ao > 
small enough. □ 

Step 3: convergence of Y^. The goal of this step is to show that the sequence Yn constructed in 
the previous step strongly converges in L°°(0, oo; £2) n L^^(0, go; endowed with the norm: 

II ■ IIl°°(0,oo;£2) + ||t'^^^ • ||l°°(0,oo;£8) + II* ^ ^- IIl°° (0,oo) ) 

and that the limit V solves the integral formulation ()4.2p of the Navier-Stokes equations p.lip - ljl.lSp . 
The main idea here comes from [15]: let us define 







Wn+l{t) :=Yn+l{t)-Yn{t) 

S{t - S)¥dw 1^, ((4 - Yn){s) {Yn - Yn-l){s) + (4 " 4-1 + Y^.i - y„)(s) ® y„_l(s)) ds. 

Again, we construct a sequence a„ such that for all n, 

a.„ > max<^ suptt || W„(t)||£8, sup || W„(t)||£2, supt^\iw„{t)\ \ ■ 

L i>0 t>0 t>0 ) 

Indeed, we have: 



3 



tmWn+limcs < K4(8,4)ii^ (t-s)"i(||y„(s)||£8 + ||y„_i(s)||£8)||W^„(s)||^8ds 

+ [\t - srHK{s)\\\Wn(,s)\\c^ + \en{s) - in-imyn-l(.s)\\c^)ds 

Jo 



< 4i^4(8,4) / (t - s) is 4^oo„ds 
Jo 



Here and in what follows, we always estimate L^-norms by interpolating the L^-norm and L^-norm 
(see (|4.8p ). In the same manner, we have 

\\Wn+imc2 < K^{2,2) [\t-s)-H\\Yn{s)\\c. + \\Yn^l{s)\\c*)\\Wn{s)\\c4ds (4.10) 

Jo 

+7^4(2,2) / (t-s)"5(|4(s)|||VF„(s)||£2 + |4(s)-4-l(s)|||yn-l(s)||£2)ds, 

Jo 
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which imphes: 

\\Wn+i{t)\\c2 < 4i^4(2,2) / {t-s)-h-^i2oands 

Jo 

< 4Cofioan. (4.11) 

Finally, we have: 

t-2\en+l{t)-£n{t)\ < i^,(4)t^(^\i-s)-t(||y„(s)||£8 + ||yn-l(s)||£8)||VF„(s)||£8ds (4.12) 

+ J\t - s)-H\in{s)\\\Wn{s)\\c4 + lUs) " 4-1 (s) 1 1| ^n-l (s) ||/:4 ) ds) 
1 /'* _3 _3 

< 4:Ki{4)t2 / {t — s) 4s 4^Qa„ds 

Jo 

< ACoHoa-n- 

Therefore, we can take a„ = (4Co/Uo)"~^ai where ai can be easily estimated thanks to Lemma 14.21 
According to Lemma [4.21 again, one can choose Aq > such that /xq < l/(4Co). With this choice, 

Edl^nll (0,00 

)) converges uniformly and there exists 

a function V £ L~(0, oo; C^) D L~ (0, oo; C^) such that 

Yn^V Strongly in L~(0, oo; C^) n L^,(0, oo; £8) , ^ iy in L^^{0, oo) . 
By construction V satisfies the decay estimates of Lemma 14.21 

snp{tl\\Vit)\\cs} < fio, sup{\\Vit)\\c2} < fio, sup{t^|^y(t)|} < /xq. (4.13) 

t>0 t>0 t>0 

The last point to check is that V indeed is a solution of the integral equation ()4.2p . i.e. we have to 
check that ICYn — t- ICV. This computation is exactly the previous one: 

)CV{t) - KYnit) = j^S{t- s)Pdiv Ir, {{Iv -V)(^{V- y„) + (^y -^n + Yn-V)® y„) [s) ds. 

3 1 

Doing as in ()4.10p and using that sup^^g ll^~^n(Oll£'* > sup^^o ll^(0~^n(i)ll/;2 and sup^^gji 2 \^y{t) — 
£„(t)|} tend to zero as n — )• 00, one easily shows that K,V — KY^ converges to in L°°(0, 00; £^). This 
shows that the limit V of the sequence Y^ solves the integral formulation ()4.2p of the Navier-Stokes 
equations (fLTB - lfrTHI) . 

Step 4: The limit V is the unique weak solution of (frTT]) - (frT8]l when Vq G H'^{W^). In the 
previous steps, we have constructed a solution to the integral formulation ()4.2p of the Navier-Stokes 
equations ()l.lip - ()1.18p verifying the — L'^ decay estimates ()4.13p . The last point that we have to 
check is that this solution V is the unique solution from the well-posedness theory of [22]. In |22j . 
uniqueness is obtained in the framework V G L°°(0,r;£2) n L'^{{),T;H'^{W^)). 

Of course, our solution satisfies by construction V G L°°(0,T;£^), and thus we only have to check 
that V G L2(0,T;i?i(M2)). 

We focus on the case of initial data Vq G H^{^^) (i.e., Vq G V{A^/'^), see ()2.ip ). In that case, we 
prove that the solution V constructed above is the unique solution in L°°(0, T; £^) nL^(0, T; H^{^)). 
The main issue is to show that the sequence ||Vyn||L°o(o,T;L2(j"o)) uniformly bounded in n for any 
arbitrary T > fixed. For T fixed, it is proved in [22l Cor. 4.3] that S{t)Vo belongs to C([0, T]-H^{R^)) 
when Vb G i7^(M^), which implies that there exists Jo > such that 

l|Vyo||L°o(0,T;L2(jrQ)) < Jq. 

Next, we construct by induction a sequence J„ such that for all n G N 

l|Vyn||L°°{0,T;L2(Jfo)) - "^n- 



36 S. ERVEDOZA, M. HILLAIRET & C. LACAVE 

Using (fLMl) with p = 2, g = 8/5 and p = ^ = 2, for all t G [0, T] 



l|Vy„+i(t)||L2(^,) < Jo + Cs/5K2{2,8/5) (i - s)-2+"8 |||y„(s)||Vy„(s)|||i8/5(^„) 



+^2(2, 2) [\t - s)-^|4(5)|||Vy„(s)||i2(^„) ds 

JO 

< Jo + Q/5^2(2,8/5) / (i-s)-trn(s)||£8||Vy„(5)||i2(^„)ds 







+i^2(2,2) / {t - s)"\en{s)\\\Vyn{s)h2iTo)ds 
Jo 

r-t 

< Jo + C8/5/^2(2,8/5) / {t- s)-ls-ifioJnds 

Jo 

+^2(2, 2) / {t- s)-h-^fioJnds 
Jo 

< Jo + CofloJn ='■ Jn+1, 



where Cg/g := 11^8/5 lli?,(L8/5(R2)^£8/5) and Co := (^8/57^:2(2, 8/5) + ^2(2, 2)) - r) It 2 dr. 
Taking Aq > small enough so that C0//0 < 1/2, there holds: 



— Joy^(CofJ,o) < -^0- ;s; < 2Jo. 

Hence we have, for all n G N, 

l|Vy„||Loo(o,T;L2{JP'o)) - ^-^0, 

which implies that verifies the same estimate. Inside -6(0, 1), we have 



Yn ■ X dx 

5(0,1) 



< CWYr 



n\\C^ 



\VYn\ = |wy„| = 

which is uniformly bounded in time and n. 

Note that this is not enoug h to conclude that Yn G L~([0, T]; i?^(E2)) for all n and one should be 
careful that the boundary conditions are compatible on OBq. In order to do that, for all n G N, we 
introduce, for all e G (0, 1) and t > 0, 

rt{l~e) 



= S{t)Vo + / ' S{t- s)FF{Ynis)) ds. 
Jo 



Of course, arguing as above, Y^^^ satisfy exactly the same estimate as Y^+i, uniformly with respect 
to e > and n G N: 

TT 

I|V5^^+i||loo(o,T;L2(R2)) = ^ll'^y„%JlL°°(0,T) + l|Vy^+i||Loo(o,T;L2(Bo)) < C. (4-14) 

But, since the semigroup S{t) is analytic on for t > 0, S{t)Vo G ^{A) (see ()2.1|) ) and 

S{t- s)FF{Ynis))ds = S{te) / 5(t(l - e) - s)PF(y„(s)) ds 
Jo 

= Site) {Yn+i{t{l - e)) - S{t{l - e))Vo) 
= S{te)Yn+i{t{l-e))-S{t)Vo. 

Since for all t > 0, Yn+i{t) G this imphes that for all t > 0, l^^+i(t) belongs to P(A) for all t > 0. 
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Besides, as e — )■ 0, Y^+i converges to Yn-^i in L^^([0, oo); since 
\\Y^^,it)-Yn+iit)\\c2 

Jt(l~e) 

< t'/^(^Ki{2,4/3)fioJol^' (l-T)-i/V-V4dr + J^i(2,2VoJo^' r'^^' dry 

Hence Yn+i is the strong limit in L°°((0, T); of the sequence of functions Y^_^_^ satisfying (j4.14p and 
the fact that for ah e > and t > 0, Y^^^{t) G V{A). Therefore, Y^+i belongs to L'^{[0,T]; H^{M.'^)). 

Besides, since the bound in (|4.14|) is uniform in e > and n S N, y also belongs to L'^{{0, T); H^(M?)). 
According to [22], when the initial data Vq belongs to H^(M?), the solution V constructed in the above 
steps is the unique weak solution of (|1.11|) - (|1.18|) . 

Step 5: Sensitivity of V to the initial data. So far, given Vq £ C'^ satisfying the smallness 
condition (j4.ip . we have constructed a solution V of the integral equation (j4.2p . In this step, we show 
that the map Vq ^ V is continuous from the ball of with radius Aq to L°°((0, oo); 

Let us consider Vq' and Vq* two elements of satisfying the smallness condition (j4.ip , and Y^ and 
Y^ the corresponding sequences in (j4.3p . We set Z„ = — Y^, which satisfies by construction 



= S{t){Vo^ - Vo') + ICY^it) - ICY^it) = Zo{t) + ICY^it) - /Cy„^(t). 
Similarly as in Step 3, we are going to construct a sequence bn such that for all n, 

5„ >max|sup{t3/8||Z„(t)||£8}, sup{||Z„(t)||£2} , snp{\tkzAt)\}\ . 



t>0 t>0 t>0 



Of course, by Theorem ll.il one can take 6o proportional to HVg" — ||£2. Since 

ICY^it) - KY^{t) = £ S{t - s)Pdiv V„ ((C - Y;^){s) Y:{s) - it - Y::){s) Y^{s)) ds 

s{t - s)Pdiv v„ ((y,^ - y:){s) ® y„^.) + Y^is) (Yj^is) - y:{s)) 
+ iC - ® Y^is) + ii{s) ® {y: - d., 

arguing as in Step 3, 

tl\\JCY:{t)-}CYM\cs < i^4(8,4)tlf At-.)-i(||y„»(.)||^s + ||y„^(.)||^8)||z„(.)||^8ds 



+ (t-.)-8(|£;;(.)|||z„(.)||^4 + |C(^)-^^(s)IK(^)llz:4)d.j 

3 /■* f 3 

< 4i^4(8,4)t8 / {t - s)~is~^iiobnds 
Jo 

where /io is the constant in Lemma 14.21 And similarly as in ()4.10p . 

\\}CY^{t)-JCYM\c2 < K,{2,2) [\t-srH\\Y:{s)\y + \\Y::{s)\\c.)\\ZAs)\\c^ds 

Jo 

+K,{2,2) r(i-s)"^(i^:;(s)iiiz„,(.)ii^2 + ic(s)-^^(^)irrf(s)ii£2)d., 

Jo 

r-t 

< 4X4(2,2) / {t- syh-^nobnds 

Jo 

< ACofJ-obn- 
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Finally, we prove as in ()4.12p that: 



+ fit - .)-^(i^:;(.)iiiz„(.)ii^4 + \c{s) - t{s)\\\Y:{s)\\c.) ds 

Jo 

< 4K^(4)t5 / (t-s)-ts-i/io6nds 

Jo 

< ACofJ-obn- 

We can then chose bn+i = bo + ACofiobn and thus (recall that 4Co/io < 1 for our choice of Aq since 
Step 3) 

In particular, passing to the limit n — t- oo, we obtain 



sup WVit) - V%t)\\c2 < C\\Vo^ - Vo'\\c2. 
t>o 

Thus, our above construction yields a map Vq ^ V continuous on the ball of of radius Aq 
to L°°((0, oo); vC^), which coincides with the map Vq ^ Vw for initial data in {{^(M?), where 
denotes the weak solution of (|l.lip - (jl.l8p . Since both maps are continuous (see "The existence part" 
in the proof of Proposition 2.5 in |22l Section 6] for the continuity of the map Vq i-^ Vw from 
to L°°((0, oo); they coincide on the ball of of radius Aq- This implies that the solution V 
constructed in Step 3, as the limit of the sequence Yn, actually is the unique weak global solution of 

Step 6: Estimates on the norm of V for 2 < p < oo. The goal of this step is to show the 
following Lemma: 

Lemma 4.3. Let Aq the constant of Lemma \4.^ For all p G [2,oo), there exists a constant H{p,Xq) 
such that for all Vq G £^ satisfying (j4.ip . the solution V of p.lip - (|1.18p satisfies: 

snp{t-^--v\\V{t)\\cp} < H{p,Xo). (4.15) 

t>o 

Proof. For p < 8 we obtain (j4.15p by interpolation of the estimates of Lemma 14.21 Assume now 
p£ [8,oo): 

t-^'-p\\V{t)\\cP < K,{p,2)\\Vo\\c2+t"^--pK^{p,4) [ {t-s)--^^-p\\V{s)\\lsds 

Jo 

ft 

+t5-?i^4(p,4) / {t-s)-'i^-p\£v{s)\\\V{s)\\c*ds 
Jo 

< Ki{p,2)\\Vo\\c2+t'^-pKi{p,A) {t-sy-^+ps--^ {fil + fil) ds 

Jo 

< Ki{p,2)Xo + 2K^{p,4:)fil (l-r)~^+PT-tdT, 

Jo 

which gives the desired estimates (|4.15p and concludes the proof of Lemma 14.31 □ 

The proof of Proposition 14.11 is then completed. □ 

Remark 4.4 (Remark on the smallness condition). The smallness condition on ||Vo||£2 is not surprising, 
and such an assumption appears in a lot of articles when global well-posedness is required (see e.g. 
[14j). In dimension 2, several works ( |21l 1271 113j in the full plane and [1] in fixed exterior domains) 
show that the L^-norm tends to zero when f — )• for initial data in L^. Of course, this allows in such 
situations to get a global result for any initial data in L^ by proving only a local result for initial data 
having small L^-norm. Unfortunately, concerning the case of a moving disk in a 2D viscous fluid. 
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despite the energy estimate satisfied by the solutions of (|l.ip - ()1.8p which immediately guarantees the 
global decay of the L^-norm of the solution, it is still not clear that the L^-norm of all solutions with 
initial data in go to zero as t — )• oo. This appears to be a challenging question. 

4.2. The case of an initial data for q G (1, 2). The goal of this section is to prove Theorem 11.31 

Proof of Theorem \1.3X The proof is based on the construction done in Proposition 14.11 

Step 1. Decay estimates for p = 2 and p = 4. Let us consider again the sequence Yn constructed 
in ()4.3p . for which we already know the decay estimates of Lemma 14.21 and, by interpolation, 

sup{t'/'\\Yn{t)\\c4} < fiQ. (4.16) 

We then prove the following lemma: 

Lemma 4.5. There exists Xo{q) small enough such that for any Vq € C'^CiC'^ satisfying the smallness 
condition (j4.ip with Aq < Xo{q), there exist constants //(2, g, Vq), -fr(4, g, Vq) for which the sequence 
Yn constructed in (|4.3p satisfies 

sup{tV«-V2||y„(t)||^,} < H{2,q,Vo), sup{tV«-V4||y„(t)||^4 < H{A,q,Vo). (4.17) 

t>0 t>0 

Consequently, we have 

sup{t'/'^~'/^Vit)\\c2} < H{2,q,Vo), sup{ti/«-^/^||y(t)||^4} < ^(4,9,^0). (4.18) 

Proof. We are looking for a sequence Hn such that for all n, 

sup{tV'/-i/2||y^(t)||^,} < sup{ti/^"V4||K„(t)||^4 < Hn. 

t>0 t>0 

Of course, Theorem 11.11 implies that Hq can be taken as Hq = {Ki{2,q) -|- J^i(4, (7))||Volk9. 
For n G N, using Theorem 11.11 and Corollary 13.101 



t'/'^-'/^Yn+i{t)\\c2 < K,{2,q)\\Vo\\c.+t'/''-'/' / K,{2,2){t-s)-'/^\\Yn{s)\\% + \£n{s)\\\Yn{s)\\c2) ds 

Jo 

Using the decay estimates of Lemma 14.21 and (j4.16p , 

\\Yn{s)\\l, < {floS~'/''){s'/''-'/mn)=fioHnS-'/'^, 
\in{s)\\\Yn{s)\\c2 < {^,,S-^'^){s^'^~^'mn)=lioHnS-^'\ 

we immediately deduce that 

snY>{t^/'^~^l^\\Yn+i{t)\\c^} < i^i(2,g)||yo||o+2i^4(2,2) ( /\l - t)"!/^"^'? dr) /io-ffn- 
t>o \Jo J 

Similar computations yield 

sup{ii/'^"i/^||y„+i(t)||^4} < Ki(4,g)||Vok^ +2if4(4,2) { [\l - t)-'/^t~'/'^ dr) fi^Hn. 
One can thus take 

-f^n+i := -f^o + Co{q)fioHn, 

with 

Co{q) := 2 (^K4{2, 2) (1 - r)"^^-!/^ dr + ^4(4, 2) ^ (1 " t)-^/\~^I^ dr^ 
Choosing Ao((7) > such that the corresponding /_fo given by Lemma W?2\ satisfies 

ao(g)/io < 1/2, (4.19) 
we thus immediately obtain that for all n, Hn+i < Hq + Hn/2, yielding Hn < 2Hq for all ?i G N. □ 
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Step 2. The solution V satisfies the decay estimates ()1.4ip . The proof of this resuh fohows 
the proof of Lemma [4.31 For p G [2,4], (|1.41|) can be deduced by interpolation with (|4.18|) . For 

p G [4, oo), we write 

t-^~p\\V{t)\\cP < Ki{p,q)\\Vo\\c'^ +t-^--^K4{p,2) [\t-s)-'+-^\\V{s)\\l,ds 



+U"Ki{p,2) I \t - sr'^-v\£vis)\\\V{s)\\c2 ds 



< K^{p,q)\\Vo\\c^ +t-^--pKi{p,2) [ {t-s)-^+h-Us{H{4,q,Vo)+H{2,q,Vo)) 

Jo 

< Ki{p,q)\\Vo\\a+K^{p,2){H{A,q,Vo)+H{2,q,Vo))fio [ {1 - r)"'^-^ t~Ut. 

Jo 



Step 3. The decay estimate on iv{t)- The proof of ()1.42p is very similar to the above one and 
is based on Corollary [3TTTJ for p > 2 such that 1/p — 1/q > —1/2, 



t'/'\iv{t)\ < K,{oo,q)\\Vo\\c',+t'^'^ Ke{p){t-s)-'/'-'/P{\\V{s)\\%, + \iv{s)\\\V{s)\\cp) ds 

Jo 

< K,{oo,q)\\Vo\\c'> 

+K,{p){H{2p, q, Vo)H{2p, 2, Vo) + H{p, q, Vo)fio) [\l - r)-V2-i/p^i/p-i/2-i/9 dr. 

Jo 

Step 4. On the map q i— >■ Xo{q). We remark that, by construction q i— )■ Xo{q) is an increasing 
function. Indeed, condition (j4.19|) indicates that our proof of Theorem [T3] requires the result of Lemma 
14.21 with fiQ = no{q) > 0, where fio{q) is an increasing function of g G (1, 2]. Since the explicit formula 
(|4.7p and (j4.9p indicates that Aq i— )■ /Uq is a continuous increasing function, the map q i— )• Xo{q) is an 
increasing function of g G (1,2]. We also note here that Xo{q) — >• when g — >• 1, since Co{q) — )• oo. 
This concludes the proof of Theorem 11.31 □ 

4.3. Proximity with the linearized semi-group. In this last subsection, we compare the asymp- 
totic structure of solutions to the Navier Stokes and Stokes equations and prove Theorem 1 1.4[ 



Proof of Theorem \1.4\ Let Vq satisfy the assumptions of our proposition. 

As g I— )• Xo{q) is an increasing function, we note that ||Vo||£2 < Xo{q) for all q G [^,2] so that V{t) 
satisfies the decay estimates of ()1.4ip - ()1.42p for arbitrary q G [q,2] and p G [2,oo). 

According to estimate ()3.33p with p G [2, oo) and = 2, for all t > 0, 

\\V{t) - S{t)Vo\\cp < [ \\S{t-s)FdiY{{V{s)-iv{s))(^V{s))\\cpds 

Jo 

< K,{p,2) f\t-sr'^'l^{\\V{s)\\l, + \iv{s)\\\V{s)\\c^)ds. (4.20) 
Jo 

But, using ()1.4ip with p = 4 and q ^ [(7, 2] and with p = 2 and q, we get: 



and using (jl.42p . we obtain: 



sup{sV<?-V4||y(,)||^,}<^(4,g~,y„), sup{sV'?-V2||y(,)||^,} < ^(2, g~ y^). 

s>0 s>0 



SViv{s^''i\lv{s)\}<H,{q,V^). 

s>0 

Hence, for all s > and q G [g, 2], we have: 

\\V{s)\\l, + \lv{s)\\\V{s)\\c'2 < H{4,q,Vofs'/'-'/'^ + H{2,q,Vo)H,{q,Vo)s'/^-^/^. (4.21) 
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The case q G (4/3,2]: proof of ()1.45p . In that case, combining ()4.20p and ()4.2ip . taking q = q, 
we immediately obtain: 

SUpt-l/P"^/2+2/'?||y(t) - S{t)Vo\\cP < C{p,q,Vo) [\l - T)"l+l/PrV2-2/g^^^ 

t>o Jo 
where we used the fact that 1/2 — 2/q > —1 for q > 4/3. 
The case q G (1,4/3): proof of (jl.43p . Here, we write 



(t _ {\\V{s)\\% + \iv{s)\\\V{s)h2^^^^) ds 

t-s)-'+y^{\\V{s)\\% + \£v{s)\\\V{s)\\c2) ds+ I 

t/2 



f (^t-s)-'+'/m\V{s)\\% + \£v{s)\\\V{s)\\c2) ds+ [ {t-s)-'+'/m\V{s)\\l, + \ev{s)\\\V{s)\\c^) ds . 

Jt/2 



Using (|4.2ip with g = 2 for s G (0, 1) and g = g for s G (l,t/2) (recall t > 2), and using t - s >t/2 
for s < t/2, 

hit) < C{p,q,Vo)t-'+'/P (^j\-^/^ds + j'^\^'^-^'^d^ < C{p,q,Vo)t-^+'/P, 

where we used that 1/2 — 2/q < —1 so that s^/^"^/'? ds < oo. 
Using ([OT]) with q = 4/3 for s G (t/2,t), we obtain 

hit) < C{p, q, Vo) f [t - s)-'-^'/Ps-' ds = C{p, q, Vo)t-'-^'/P [' (1 - r)''+'/'PT-' dr. 

Jt/2 Jl/2 

The case q = 4/3: proof of (|1.44p . This case follows similarly as the previous one, except that 
estimating Ji yields 

hit) < C{p, q, Vo)t-'+'/P ds + S-' ds^ < C{p, q, Vo)t-'+'/P{l + log(t)). 

We now concentrate on the estimate (|1.46p - ()1.48p on ivit) — (^s{t)Vo- order to do that, again, we 
split the integral in two parts: 

\Mt) - %hJ < Ke{2) j'^\t - s)~^(\\V{s)\\% + \iv{s)\\\V{s)\\c2) ds 

+IUp) f\t - s)-^--^i^ {\\V{s)\\%, + ms)\\\V{s)\\c.) ds =: J^{t) + J2{t) . 

where Kg is the constant of Corollary 13 . 1 1 1 and p G (2, oo). The estimate of Ji can be done as previously 
by using (|4.2ip : 

r C{p,q,Vo)t^/^-^/'i ifgG (4/3,2], 

Ji{t)<{ C{p,q,Vo)t-Hl + log{t)) ifg = 4/3, 

[ C{p,q,Vo)t-^ if gG (1,4/3) 

For J2, remark that similarly as in ()4.2ip we can obtain for all G [q,2], s > 0, 

\\V{s)\\%, + \iv{s)\\\V{s)\\cP < {H{2p,q,Vo))h'/P~^/'^ + H{p,q^^^^ 

so that: 

J2{t)<C{p,q,Vo)t'/'-'/'i. 
This ends the proof by choosing q = q for q £ (4/3, 2] and q = 4/3 if (/ G (1,4/3]. □ 
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5. Further comments 



We list below several comments. 

Concerning optimality of Theorem When considering the decay estimates of Theorem it is 
natural to ask if the results are sharp, in particular regarding the decay of the gradient estimates when 
p > 2 and t > 1, since all other decay estimates correspond to the classical ones for the heat semigroup 
on M^. However, in our proof, the decay estimate (jl.35p differs from the one corresponding to the heat 
semigroup on M? for all the modes. Each time, this slower decay rate for t > 1 and p > 2 arises due 
the presence of the boundary. Let us point out that P. Maremonti and V. A. Solonnikov prove in [T7j 
that, when considering the Stokes equations in an exterior domain of with homogeneous Dirichlet 
boundary conditions, estimate ()3.3p . which is the counterpart of ()1.35p . is sharp. It is then likely that 
estimates of Theorem 11.11 are sharp as well. 

Straightforward extensions of theorems \1.4\ o-n-d 

• Using the density of H in n £^ and the decay estimates of Theorem ll.H one easily get, 
for all q G (1,2], for all Vq ^ Cfl ^ C? satisfying ||Vo||£2 < Ao(5/4), the unique associated solution Vif) 



Indeed, for Vb G ^^n/:^ satisfying HFolIra < Ao(5/4), by TheoremOl _ S'(t)Vb||£p goes 

to zero as t — CO (recalling that Ao(g) > Ao(5/4) if q > 5/4, see the end of Introduction). We then 
take e > and Vb G C}r\C? satisfying \\Vq - %\\ci < e. According to Theorem fTTl t^/'i-^/P\\S{t)Vo - 
S{t)Vo\\cp < Ki{p,q)u^/P~^lie. But Theorem O also implies \\mt^^t^/i~^IP\\S{t)Vo\\ = since Vq 
belongs to for some g E iX^l)- Hence 



Since e was arbitrary, this implies ()5.ip . 

• The proofs of Theorems I1.3H1.4I are only based on the — L'^ estimates for the Stokes problem 
given in Theorem 11.11 As such estimates are already known in the case of a fixed exterior domain (see 
[HlEllTT]), we claim that Theorem 11.41 holds true also in this case. Hence, the computations herein 
extend the results in [101 [T2] to the case of finite energy initial data. 

• In order to obtain the decay estimates of Theorem II. H our approach is strongly based on the 
fact that the rigid body is an homogenous disk. Indeed, in polar coordinates we decompose in Fourier 
series. A case which can be easily treated by our analysis is when the disk is non-homogenous, and 
the center of mass corresponds to the center of the disk (e.g. for a density p with radial symmetry). 
In this case, the equations (|l.lip - ()1.18p and ()1.22p - (|1.29p are the same, where: 



To our knowledge, the case of a more general shape or more general density is completely open. A 
similar problem, also open, would be to derive decay estimates in the case of two rigid disks. 

Open problems. 

• Despite Theorem [T31 a complete description of the first term in the asymptotic behavior as t — )• oo 
of the solutions V of ()l.lip - ()1.18p is still missing. Indeed, Theorems 11.21 and 11.41 cannot be combined 
since Theorem 11.21 requires the initial data to be and in that case, Theorem 11.41 onlv yields that the 
£P-norm of the difference between the solution of the complete non-linear system (|l.lip - ()1.18p and 
the linear one, given by 5(t)Vb, decays as Ct^/^~'^ , which is precisely the order of magnitude of the 
D'-noim. of the solution of the linear Stokes equation when Vq S >C^. At this level, let us emphasize 
that one of the main conceptual difficulties of this problem is that the invariant seems to be the C^- 
norm of the solution of ()l.lip - ()1.18p . despite the fact that the linear semigroup does not seem to be 
well-posed in C^. To justify this statement, we emphasize that the asymptotic given by Theorem 11.21 



to (frTT]l - (frT8]l satisfies: 



lim t^/'i-^'P\\V{t)\\cv =0. 



(5.1) 



limsupti/'?-^/P||y(t)|| < Ce. 
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does not belong to C} . Showing that the non-hnear term decreases faster than is the major 

issue which prevents us from extracting the asymptotic first order. 

If we look carefully at the proof of Theorem 11.41 we note that the difficulty comes from the fact 
that we do not manage to prove that ||S(t)Pdiv -F||£p decays faster for F ^ C with q <2 (see (|3.34|) ). 
In particular, in the case of the Navier-Stokes equations in M^, using heat kernel estimates, which are 
better than estimates (j3.34|) when q < 2 and t > 1, A. Carpio in [3] shows that the non- linear term is 
smaller than the Stokes solution for large time. But as we have noted above, the restriction in (|3.34|) 
seems to be unavoidable. 

Nevertheless, other methods relying on the use of suitable scaling invariance and similarity variables 
have been used for providing leading order terms in [3l [El [11]. To keep the unity of this paper we 
postpone these approaches to a future work. 

• Another open problem is to remove the smallness condition in Theorem 11.41 it is done for 
Navier-Stokes in the full plane [211|27] and in fixed exterior domains [1]. Indeed, such result would be 
expected in view of the energy dissipation law (jl.l9p which indicates the decay of the i2^-norm of the 
solutions of ([rTD-KLSD. 
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Appendix A. Proof of Proposition 12.21 

Assume Vq £ Ci C^{M?) Hille-Yosida's theorem implies there exists a unique solution V G 
C([0,oo);£2) to ([L22])-(IL29])- Furthermore, the unknowns {iy , cjy) and the pressure p that are con- 
structed starting from V have, with V, the following regularity (see \22\ Corollary 4.3]): 

V e C([0,oo); [H\To)]^) n L\{0,<x); [H\F^)f) , Vp G L\{Q,<^)-L\Fo)) , 

iv£H^{0,oo), ujy £ H^{0,oo). 

We note also that further smoothing properties of the semigroup (see [21 Theorem 7.7]) imply 

{v,p) G [C°°((0,oo) X Jo)]', VpGC((0,oo);L2(Jo)). 

Consequently, we introduce the decomposition (IF, <I>, Vr) of V in spherical harmonics and a 
corresponding decomposition of the pressure p: 

1 /■^'^ 

pi{t,r) = - p{t,r,e) cos 9 d9 , (A.l) 



vr 



1 







vr 



2-7T 



qi{t,r) = - / p{t,r,9)sm9d9 , (A.2) 







PR{t,r) = p{t, r, 9) — pi{t,r) cos 9 — qi{t,r) sin 9 . (A. 3) 

Note that, like vr, the remainder term pn satisfies: 



/ pB{t,r,9)cos{9)d9 = / pR{t,r,9) sm{9) d9 = 0. (A.4) 
Jo Jo 

Applying the continuity of the spherical-harmonic decomposition together with the continuity of V 

yields: 

W G C([0,oo);L2((0,oo),rdr)) (a^^,*/r) G C([0, oo); ^^((0, oo), rdr)) 
Vr G C([0,oo);L2(Jo)) (5,$,$/r) G C([0, oo); ^^((o, oo), rdr)). 
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together with 

drPi £ C((0,oo);L2((i,oo),rdr)), drQi G C((0, oo); L2((1, oo), rdr)). 

Referring to the formulas (j2.8|) - (|2.1U|) and ()A.ip - ()A.2p we also obtain at once the smoothness of 
{WjVr,^,"^), of ipi,qi) and of (VrjPji). It now remains to compute the systems satisfied by these 
unknowns. 

We recall that the spherical-harmonic decomposition of V reads V = VrCr + Vge^ 

Vr = —sm9 cos 6* + Vr • e,. , Vg = W mm(l,r) + dr'i/ cos9 + dr^ sin9 + Vr ■ eg (A. 5) 

r r 

and the velocity- field on the disk is given as follows in radial coordinates: 

{iv + ujvx~^)r = (-v,i cos{9) + lv,2 siu{9) , {iv + wyx"'")^ = ujr — iv,i sm{9) + iv,2 cos{9). 

Identifying V and the velocity-field of the disk on dBo ( i.e. for r = 1), we obtain the following 
boundary conditions: 

w{t,l) = uj{t), yt>o, 

i;{t,l) = iv,2it), dMt,l) = iv,2{t), yt>0, 

ipi{t,l) = -iv,i{t), drMtA) = -^vAt), Vt>0. 

VR{t,x) = 0, \fxedBo, \/t>0. 

In the fluid domain, we remark that, introducing x such that drX = w and x(0) = 0) the spherical- 
harmonic decomposition reads: 

V = V-^ {x + ip cos{9) + if sm{9)) + vr 

so that: 

dtv — uAv = 'V'^[dt — i^A] (x + 'ipcos{9) + ipsm{9)) + dtVR — i^Avr , 
where, in polar coordinates: 

We also recall that, the gradient operator reads: 

Vq = Orqer H eg . 

r 

Finally, we remark that orthogonality conditions such as (|2.2p or ()A.4p transmit to time and space 
derivative. Hence, replacing ip and p by their values in the two last formulas, identifying then the 
different frequencies: constant, cos9, sm9, and remainders, we get: 

dtw — V ( drrW + — ^ ) = , for (t, r) G (0, oo) x (1, oo); 
\ r ) 

drrip + ^ j ^ -fdrqi , for (t, r) G (0,oo) X (l,oo); 

dtdrtp — vdr (drrip + — ) = -— , for it, r) G (0, oo) X (1, oo); 

dtip - V (^drrf + - = rdrPi , for (t, r) G (0, oo) x (1, oo); (A.6) 
dtdr^ - udr (^drrV + ^ - = ^ , for (t, r) G (0, oo) X (1, oo); (A.7) 
dtVR - uAvR + VpR = 0, t > 0, X G Jo • (A.8) 
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To end up the proof of Proposition 12.21 we write now ()1.26p - ()1.27p . First we recall that, on Bq there 
holds n = —Br (the normal n is here computed outward the fluid domain) and r = 1, so that: 

-2D{v)n = drv + Vvr - [Ver]^v , 

= {2drVr)er + {drVg + dgVr - Ve)ee , 
— En = {—p + 2l'drVr)er + u{drVg + dgVr — Vg)eg. 

Hence, computing for instance = ^'v ' ^2 we have: 

/'27r r2n 
"^^'^,2 = / {2vdrVT — p) sin 6d9 + U I {dr-Vg + dgVr — Vg) COS 6 d9 , 

Jo Jo 

= / {2l' drVr + — p)svil 6 AO + U j {drVg + dgVr) COS 9 d9 . 

Jo Jo 
In these last integrals, we then compute drVr and drVg with respect to w, ip and tp, and vji thanks to 
(|A.5p . Recalling the orthogonality conditions (j2.2p and (|A.4p . we get: 



V,2 



vr 



vr 



2dr 



r 



+ drrip + 1p - drip ] - Ql 



V OrrW H 9 -^1 



The computations of Jio'y and vn^y-^ are similar. 
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